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ABSTRACT 

We perform a comprehensive estimate of the frequency of galaxy mergers and their impact on star 
formation over z - 0.24-0.80 (lookback time Tb ~ 3-7 Gyr) using - 3600 (M > 1 x 10^ Mq) galaxies 
with GEMS HST, COMBO-17, and Spitzer data. Our results are: (1) Among - 790 high mass 
(M > 2.5 X 10^° ^'^o) galaxies, the visually-based merger fraction over z ~ 0.24-0.80, ranges from 9% 
± 5% to 8% ± 2%. Lower limits on the major merger and minor merger fraction over this interval range 
from 1.1% to 3.5% , and 3.6% to 7.5%, respectively. This is the first, albeit approximate, empirical 
estimate of the frequency of minor mergers over the last 7 Gyr. Assuming a visibility timescale of 
^ 0.5 Gyr, it follows that over Tb ~ 3-7 Gyr, ~ 68% of high mass systems have undergone a merger of 
mass ratio > 1/10, with ^ 16%, 45%, and 7% of these corresponding respectively to major, minor, and 
ambiguous 'major or minor' mergers. The average merger rate is a few xlO"** galaxies Gyr~^ Mpc~'^. 
Among ^ 2840 blue cloud galaxies of mass M > 1.0 x 10^ M©, similar results hold. (2) We compare the 
empirical merger fraction and merger rate for high mass galaxies to three ACDM-based models: halo 
occupation distribution models, semi-analytic models, and hydrodynamic SPH simulations. We find 
qualitative agreement between observations and models such that the (major -l-minor) merger fraction 
or rate from different models bracket the observations, and show a factor of five dispersion. Near- future 
improvements can now start to rule out certain merger scenarios. (3) Among ~ 3698 M > 1.0 x 10^ Mq 
galaxies, we find that the mean SFR of visibly merging systems is only modestly enhanced compared to 
non-interacting galaxies over z ^ 0.24-0.80. Visibly merging systems only account for a small fraction 
(< 30%) of the cosmic SFR density over Tb ~ 3-7 Gyr. This complements the results of Wolf et al. 
(2005) over a shorter time interval of Tb ~ 6.2-6.8 Gyr, and suggests that the behavior of the cosmic 
SFR density over the last 7 Gyr is predominantly shaped by non-interacting galaxies. 

Subject headings: galaxies: fundamental parameters — galaxies: structure — galaxies: kinematics and 
dynamics — galaxies: evolution 



1. INTRODUCTION 

Hierarchical A cold dark matter (ACDM) models pro- 
vide a successful paradigm for the growth of dark matter 
on large scales. The evolution of galaxies within ACDM 
cosmogonies depends on the baryonic merger history, the 



star formation (SF) history, the nature and level of feed- 
back from supernovae and AGN, the redistribution of an- 
gular momentum via bars or mergers, and other aspects of 
the baryonic physics. Empirical constraints on the fate of 
the baryonic component, in particular their merger and SF 
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history, are key for developing a coherent picture of galaxy 
evolution and testing galaxy evolution models (e.g., Kauff- 
mann et al. 1993; Somerville & Primack 1999; Navarro 
& Steinmetz 2000; Murali et al. 2002; Governato et al. 
2004; Springel et al. 2005a,b; Mailer et al. 2006). Such 
constraints can help to resolve several major areas of dis- 
cord between observations and ACDM models of galaxy 
evolution, such as the angular momentum crisis, the prob- 
lem of bulgeless and low bulge-to-total {B /T) ratio galax- 
ies (Navarro & Benz 1991; D'Onghia & Burkert & 2004; 
Kautsch et al. 2006; Barazza, Jogee, & Marinova 2008; 
Weinzirl et al. 2009), and the substructure problem. 

The merger history of galaxies impacts the mass assem- 
bly, star formation history, AGN activity and structural 
evolution of galaxies. Yet, the merger rate has proved 
hard to robustly measure for a variety of reasons. Initially, 
small samples hindered the first efforts to measure merger 
rates (Patton et al. 2000; Le Fevre et al. 2000; Conselice 
et al. 2003). Later studies drew from larger samples, and 
have used a variety of methods to characterize the inter- 
action history of galaxies at 2 < 1. Studies based on close 
pairs report a major merger fraction of 2% to 10% over 
z~ 0.2 to 1.2 (Lin 2004; Kartaltepe et al. 2007; de Ravel 
et al. 2008) and ^ 5% for massive galaxies at z~ 0.4 to 
0.8 (Beh et al. 2006). 

Among these studies, the estimated fraction of galaxies, 
which have experienced a major merger since z 1.2 dif- 
fer. Lin et al. (2004) find that ^ 9% of the massive early 
type galaxies experienced a major merger since z ~ 1.2, 
while Bell et al. (2006) report that ~ 50% of all galaxies 
with present-day stellar masses above 5 x 10^*^ Mq have 
undergone a major merger since z ~ 0.8. This difference 
of a factor of five between the two studies is addressed by 
Bell et al. (2006; see their footnote 9) and is traced pri- 
marily to differences in the following factors: the dataset 
used, the redshift integration method, the way in which 
the authors handle the relative number density of the par- 
ent population from which pairs are drawn and that of the 
remnant population, the assumed visibility timescale, and 
the fraction of pairs estimated to be true gravitationally 
bound pairs. 

Studies based on Gini-M20 coefficients report a fairly 
constant fraction (^ 7 ± 0.2 %) of disturbed galaxies over 
z^ 0.2 to 1.2 among bright galaxies (Lotz et al. 2008) in 
the AEGIS survey. Similar trends are found from early 
results based on visual classification and asymmetry pa- 
rameters of high mass galaxies (e.g., Jogee et al. 2008). 
The study by Cassata et al. (2005) based on both pairs 
and asymmetries report a mild increase in the merger rate 
with redshift up to z~ 1, with the caveat of a small sample 
size. The merger rate/fraction at z > 1 remains highly un- 
certain, owing to relatively modest volumes and bandpass 
shifting effects, but there is general trend towards higher 
merger fractions at higher redshifts (e.g., Conselice et al. 
2003; Cassata et al. 2005). 

Studies to date have brought important insights but face 
several limitations. In the case of studies based on close 
(separation ~ 5 to 30 kpc) pairs, the translation of the pair 
frequency into a merger rate is somewhat uncertain due to 
several factors. The uncertainties in the spectrophotomet- 
ric redshifts for some of the galaxies in pairs can cause us 

We assume a flat cosmology with Om,o = 1- ^^A = 0-3 and Ho 



■ overestimate or underestimate the true close pair frac- 
tion, with the latter effect being more likely. Corrections 
for this effect are uncertain and depend on the shape of 
the spcctrophotometric redshift errors (e.g., see Bell et al. 
2006 for discussion) . Secondly, even pairs with members at 
the same redshift may not be gravitationally bound, and 
may therefore not evolve into a merger in the future: this 
effect causes the close pair fraction to be upper limits for 
the merger fraction. Thirdly, gravitationally bound pairs 
captured by this method sample different phases of an in- 
teraction depending on the separation, and any merger 
rate inferred depends on the separation, orbital eccentric- 
ity, and orbital geometry. 

In the case of studies, which use automated parame- 
ters, such as CAS asymmetry A and dumpiness param- 
eters (Conselice et al 2000; Conselice 2003) or Gini-M20 
coefficients (e.g., Lotz et al. 2004) to identify merging 
galaxies can fail to pick stages of both major and minor 
mergers where distortions do not dominate the total light 
contribution (§[3T4]). Comparison with simulations suggest 
that the CAS criterion {A > 0.35 and A > S; Conselice 
2003) capture major mergers about 1/3 of the time, while 
the eye is sensitive to major merger features over twice as 
long (e.g., Consehce 2006; i; 13. 4p . To complicate matters, 
automated asymmetry parameters can also capture non- 
interacting galaxies hosting small-scale asymmetries that 
are produced by stochastic star formation ('? l4.2p . Thus, it 
is important to use several methods to assess the robust- 
ness of results and identify the systematics. 

In this paper, we present a complementary study of the 
frequency of mergers and their impact on the SF activity 
of gala xies over z ^ 0.24-0.80 (lookback times Tb of 3-7 
GyilE]) using HST ACS, COMBO-17, and Spitzer 24 /im 
data of 3600 galaxies in the GEMS survey. The outline 
of the paper is given below and describes how this study 
complements existing work: 

1. We use a large sample of ~ 3600 (M > 1 x 10^ 
Mq) galaxies to get robust number statistics for the 
merger fraction among ~ 790 high mass (M > 2.5 x 
10^° Mq) galaxies and - 2840 M > 1 x 10^ Mq blue 
cloud galaxies (§[1 Tabled Tabled). 

2. Two independent methods are used to identify 
merging galaxies: a physically-driven visual clas- 
sification system complemented with spcctrophoto- 
metric redshifts and stellar masses 13.21 to i) I3.3p . 
as well as automated CAS asymmetry and dumpi- 
ness parameters (ij l3.4p . This allows one of the most 
extensive comparisons to date between CAS-based 
and visual classification results I4.2p . 

3. We design the visual classification system in a 
way that allows merger fractions and rates from 
observations and theoretical models to be read- 
ily compared. We classify as mergers those sys- 
tems that show evidence of having experienced a 
merger of mass ratio > 1/10 within the last visibil- 
ity timescale. Throughout this paper, we use the 
standard definition whereby major mergers are de- 
fined to have stellar mass ratio (1/4 < M1/M2 < 
1/1), while minor mergers have (1/10 < M1/M2 < 
1/4). We set lower limits on the major and mi- 
nor merger fraction 14. ip . To our knowledge, this 

km s~^ Mpc^^ throughout. 
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is the first, albeit approximate, empirical estimate 
of the frequency of minor mergers over the last 7 
Gyr. While many earlier studies focused on major 
mergers, it is important to constrain minor merg- 
ers as well, since they dominate the merger rates in 
ACDM models, and play an important role in build- 
ing the bulges of massive galaxies (e.g., Weinzirl et 
al. 2009). 

4. We compare the empirical merger fraction and rate 
to a suite of ACDM-based simulations of galaxy 
evolution, including halo occupation distribution 
models, semi-analytic models, and hydrodynamic 
SPH simulations 14. 5p . To our knowledge, these 
extensive comparisons have not been attempted to 
date, and are long overdue. 

5. The idea that galaxy interactions generally enhance 
the star formation rate (SFR) of galaxies is well es- 
tablished from observations (e.g., Larson & Tins- 
ley 1978; Joseph & Wright 1985; Kennicutt et al. 
1987; Barton et al 2003) and simulations (e.g., Ne- 
groponte & White 1983; Hernquist 1989; Barnes 
& Hernquist 1991, 1996; Mihos & Hernquist 1994, 
1996; Springel, Di Matteo & Hernquist 2005b). 
However, simulations cannot uniquely predict the 
factor by which galaxy mergers enhance the SF ac- 
tivity of galaxies over the last 7 Gyr, since both the 
SFR and properties of the remnants in simulations 
are highly sensitive to the stellar feedback model, 
the bulge-to-disk (B/D) ratio, the gas mass frac- 
tions, and orbital geometry (e.g.. Cox et al 2006; 
di Matteo et al. 2007). This motivates us in § ITBl 
to empirically investigate the impact of interactions 
on the average UV-based and UV-l-IR-based SFR of 
high mass (M > 2.5 x 10^° Mq) and intermediate- 
to-high mass {M > 1 x 10^ -^0) galaxies over 
z - 0.24-0.80. 

6. The SF properties of merging and non- interacting 
galaxies since z < 1 is of great astrophysical inter- 
est, given that the cosmic SFR density is claimed 
to decline by a factor of 4 to 10 since z ~ 1 (e.g., 
Lilly et al. 1996; Ellis et al 1996; Flores et al. 1999; 
Haarsma et al. 2000; Hopkins 2004; Perez- Gonzalez 
et al. 2005; Le Floc'h et al. 2005). In § mil we 
set quantitative limits on the contribution of visi- 
bly merging systems to the UV-based and UV+IR- 
based SFR density over z ~ 0.24-0.80. Our study 
covers a 4 Gyr interval (Tb ~ 3-7 Gyr or 2; ~ 0.24- 
0.80) and extends the earlier studies carried out 
over a smaller 0.6 Gyr interval (Tb ~ 6.2-6.8 Gyr 
or z ~ 0.65-0.75) by Wolf et al. (2005) and Bell 
et al. (2005) on the UV and IR luminosity density, 
respectively. Our study also complements IR-based 
studies by Lotz et al. (2008), Hammer et al. (2005; 
195 galaxies at z > 0.4 ) and Melbourne et al. (2005; 
~ 800 galaxies) in terms of better number statistics 
and the use of both UV-based and IR-based SFR 
indicators. 

2. DATASET AND SAMPLE SELECTION 

This study uses data from the Galaxy Evolution from 
Morphology and SEDS (GEMS; Rix et al. 2004) survey. 

We adopt a Chabrier (2003) IMF when exploring the contribution 



which provides high resolution Hubble Space Telescope 
(HST) Advanced Camera for Surveys (ACS) images in 
the F606W and F850LP fihers over an 800 arcmin^ (~ 
28' X 28') field centered on the Chandra Deep Field-South 
(CDF-S). Accurate spectrophotometric redshifts [(5z/(l -I- 
z) ~ 0.02 down to i?vega — 24] and spectral energy distri- 
butions, based on 5 broad bands {U BVRI) and 12 medium 
band filters, are available from the COMBO-17 project 
(Wolf et al. 2004). The ACS data reach a limiting 5 cr 
depth for point sources of 28.3 and 27.1 AB mag in F606W 
and F850LP, respectively (Rix et al. 2004; CaldweU et al. 
2008). The effective point spread function (PSF) in a sin- 
gle F606W image is ~0'.'07, corresponding to 260 pc at 
z ~ 0.24 and 520 pc at z - 0.80. The PSF of combined 
drizzle images is ~0'.'l. In addition to HST ACS imag- 
ing, the GEMS field has panchromatic coverage which in- 
cludes Spitzer (Rieke et al. 2004; Papovich et al. 2004) and 
Chandra data (Alexander et al. 2003; Lehmer et al. 2005). 

We use stellar masses from Borch et al. (2006). We refer 
the reader to the latter publication for a detailed descrip- 
tion and provide a summary here. Using the 17-passband 
photometry from COMBO-17, objects were classified as 
main sequence, stars, white dwarfs, galaxies, and quasars 
using color indices and their photometric redshifts were 
estimated using simple dust-reddened single-burst SED 
templates (Wolf et al. 2004). For galaxies and quasars, 
the joint probability of a given redshift and a given rest- 
frame SED is derived and this procedure provides a min- 
imum error variance estimation of both the redshift and 
the SED template. Once the redshfit has been estimated, 
the SEDs in 17 bands were fitted with a new set of tem- 
plate SEDs with more plausible SF histories in order to 
derive a stellar M/L (Borch et al. 2006). The library of 
SEDs is built using the PEGASE stellar population syn- 
thesis model (see Fioc & Rocca-Volmerange 1997 for an 
earlier version of the model) and the underlying SF his- 
tories are parameterized by the three-component model, 
with a Kroupa (Kroupa et al. 1993) initial mass function 
(IMF) adopted in the mass regime 0.1-120 Mq. We note 
that the stellar masses are consistent within 10% with the 
masses that would be estimated using a different Kroupa 
(2001) or Chabrier (2003) IMF0. The reddest templates 
have smoothly- varying exponentially-declining star forma- 
tion episodes, intermediate templates have a contribution 
from a low-level constant level of star formation, while the 
bluer templates have a recent burst of star formation su- 
perimposed. Random stellar mass errors are < 0.3 dex on 
a galaxy-by-galaxy basis, and systematic errors in the stel- 
lar masses were argued to be at the 0.1 dex level (Borch 
et al. 2006). Beh & de Jong (2001) argued that galax- 
ies with large bursts of recent star formation could drive 
down stellar M/L values by up to 0.5 dex at a given color, 
but we note that the Borch et al. (2006) templates do in- 
clude bursts explicitly, thus compensating for the worst 
of the uncertainties introduced by bursting star formation 
histories. 

We present the results based on the visual classifica- 
tion and CAS (Conselice 2003) parameters (§[2|) of GEMS 
F606W, rather than F850LP images, for the following rea- 
sons. The F606W images are ~ 1.2 AB magnitude deeper 
than the GEMS F850LP images and allow more reliable 

mergers to the cosmic SFR (ij 14.61 and ij 14.71 1. 
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characterization of morphological features in the presence 
of cosmological surface brightness dimming at the rate of 
(l+z)"^ (e.g-, Barden et al. 2007). Furthermore, the low 
signal to noise in the F850LP images leads to large error 
bars on the asymmetry A and dumpiness S parameters 
generated by the CAS code, effectively making it imprac- 
tical to use these values in CAS merger diagnostics (i) l3.4p . 
When using the F606W images, we only include results 
over the redshift range z 0.24-0.80 in order to ensure 
that the rest-frame wavelength Arost stays primarily in the 
optical band and does not shift well below the 4000 A 
break. In the fourth redshift bin (z ^ 0.6 to 0.8) Aicst 
shifts to the violet/near-UV (3700 A to 3290 A), but as 
we show in § 14.11 this does not significantly impact the 
results. We discard the last redshift bin at z > 0.8 where 
Arest shifts into the far-UV. These steps lead to a sam- 
ple of 4766 galaxies selected down to -Rvoga < 24, over 
z ^ 0.24-0.80 (Tback ^ 3-7 Gyr). 

In this paper, we present results for two samples of as- 
trophysical interest, which are derived by applying stellar 
mass cuts to the above sample of ^ 4766 galaxies The first 
sample (henceforth sample SI) focuses on galaxies with 
high stellar mass (M > 2.5 x 10^° Mq; Table[T]). The sam- 
ple size is originally 804 galaxies, out of which 798 (99.2%) 
could be visually classified. For this stellar mass range, the 
red sequence and blue cloud galaxies are both complete 
out to the highest redshift bin z ~ 0.62-0.80 for our sam- 
ple, and we have theoretical predictions for comparison 
(see § 14. 5p from semi-analytical models (e.g, Somerville et 
al. 2008; Hopkins et al. 2007; Khochfar & Burkert(2005); 
Bower et al. 2006), TV-body (D'Onghia et al. 2008), and 
hydrodynamical SPH simulations (e.g.. Mailer et al. 2006). 
Note that the survey has few galaxies above 10^^ Mq 
(Fig. [J), and hence the high mass (M > 2.5 x 10^" Mq) 
sample primarily involves galaxies in the range 2.5 x 10^*^ 
to 10" Mq. 

We also present selected results for the sample S2 of 
~ 3698 galaxies with M > 1 x 10^ Mq and visual classes. 
Although this sample includes the ~ 790 high mass galax- 
ies in sample SI, it is dominated by systems of intermedi- 
ate mass (10^ < M/Mq < 2.5 x 10^°). For the mass range 
M > 1 X 10^ Mq, the blue cloud is complete in our sample 
out to z ^ 0.80, while the red sequence is incomplete in the 
higher redshift bins. Where appropriate, we will therefore 
present results for the blue cloud sample only (e.g., lower 
part of Table ^ . The rest-frame U — V color is plotted 
versus the stellar mass for the sample S2 in Figure [TJ The 
redshift interval is divided into four 1 Gyr bins. The diag- 
onal line marks the separation of the red sequence and the 
blue cloud galaxies (BCG) at the average redshift Zavo of 
the bin. We use the definition in Borch et al. (2006) and 
Bell et al. (2004) for CDF-S: 

([/ - V"),.est > 0.227 log(M/MQ) - 1.26 - 0.352z (1) 

The vertical lines on Figure [T] marks the mass complete- 
ness limit (Borch et al. 2006) for the red sequence galaxies. 
The blue cloud galaxies are complete well below this mass 
(Borch et al. 2006). 

3. METHODOLOGY: IDENTIFYING MERGERS AND 
NON-INTERACTING GALAXIES 



3.1. Overview of the Methodology 

Galaxy mergers with a mass ratio M1/M2 > 1/10 can 
have a significant impact on galaxy evolution. Accord- 
ing to simulations, major mergers (defined as those with 
mass ratio 1/4 < M1/M2 < 1/1) typically destroy stel- 
lar disks, transforming them via violent relaxation, into 
systems with a steep or r^^^ de Vaucouleurs-type stel- 
lar profile, such as ellipticals (e.g., Negroponte & White 
1983; Barnes & Hernquist 1991; Mihos & Hernquist 1996; 
Struck 1997; Naab & Burkert 2001; but sec Robertson et 
al. 2004). Simulations further suggest that recent major 
mergers at z < 1 are typically associated with arcs, shells, 
ripples, tidal tails, large tidal debris, highly asymmetric 
light distributions, double nuclei inside a common body, 
galaxies linked via tidal bridges of light, and galaxies en- 
closed within the same distorted envelope of light. 

Minor mergers (defined as those with 1/10 < M1/M2 < 
1/4) of two spirals will not destroy the disk of the larger 
spiral (e.g., Hernquist & Mihos 1995; Smith et al. 1997; 
Jogee et al. 1999). Typically, the smaller companion sinks 
via dynamical friction, may excite warps, bars, spirals, and 
other non-axisymmetric perturbations, and leads to verti- 
cal heating, arcs, shells, ripples, tidal tails, tidal debris, 
warps, offset rings, highly asymmetric light distributions, 
etc (e.g., Quinn et al. 1993; Hernquist & Mihos 1995; Mi- 
hos et al. 1995; Quinn, Hernquist, & FuUagar 1993; Smith 
et al. 1994; Jogee 1999; Jogee et al. 1998, 1999; review by 
Jogee 2006 and references therein). 

One goal of this paper is to identify systems whose mor- 
phology and other properties suggest they they have re- 
cently experienced a merger of mass ratio M1/M2 > 1/10. 
We employ two methods: a physically-driven visual clas- 
sification system (ij 13.21 to ij 13. 3p complemented with spec- 
trophotometric redshifts and stellar masses, and a method 
based on quantitative asymmetry (A), and dumpiness (5) 
parameters (§ 13. 4p derived using the CAS code (Conselice 
2003). While many studies use only automated methods 
or visual classification, we choose to use both methods in 
order to better assess the systematics, and to test the ro- 
bustness of our results. 

3.2. Visual Classification of Mergers 

The visual classification system we adopted for identify- 
ing mergers was aimed at setting a procedural framework 
that allows merger fractions and rates from observations 
and the theoretical models (outlined in § I4.5[) to be defined 
in similar ways and to be readily compared. 

Specifically the theoretical models described in § 14.51 
track systems which have a stellar mass M^, > Mcut and 
have experienced a merger of a certain mass ratio M1/M2 
within the last visibility timescale tvis , between times (tobs 
- tvis) and tohs- Here fobs is the time corresponding to the 
observed redshift z; tvis is the timescale over which mor- 
phological distortions remain visible after a merger and we 
adopt a nominal value of 0.5 Gyr (see § 14. 3p : Mcut is the 
cutoff mass for the merger, which is 2.5 x 10^° Mq for the 
high mass sample SI and 1 x 10^ Mq for the intermediate 
mass sample S2 (§[2|). In the models, major mergers are as- 
sociated with a mass ratio of 1/4 < M1/M2 < 1/1, minor 
mergers with 1/10 < M1/M2 < 1/4, and major-|-minor 
mergers with M1/M2 > 1/10. 

Analogous to the theoretical models, the goal of the 
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visual classification system is to identify systems with 
> Mcut, which show evidence of having experienced 
a merger of mass ratio > 1/10 within the last visibility 
timescale tvis- In this paper, we refer to such systems ei- 
ther as interacting galaxies or as mergers. These mergers 
are denoted as having visual type 'Int' in the figures of 
this paper and examples are shown in Figure [H In prac- 
tice, they consist of three types of mergers, which we de- 
note as Int-1, Int-2a, and Int-2b, and define in 13.2. II and 
§ 13.2.21 These three types of systems encompass young to 
advanced mergers, and are identified/handled in different 
ways, as described below. 

3.2.1. Mergers oj type Int-1 

Mergers of type Int-1 primarily represent advanced 
mergers of mass Af* > Mcut, which appear as a single 
system in ACS images, and are likely associated with a 
merger of mass ratio > 1/10 that occurred betwen times 
(iobs - ^vis) and tobs- Systems of type Int-1 are identified 
empirically based on the following criteria: 

1. They have M, > Mcut and show morphological dis- 
tortions, which are similar to those seen in simula- 
tions of mergers of mass ratio > 1/10, and remain 
visible over the visibility timescale ivis- The distor- 
tions include arcs, shells, ripples, tidal tails, large 
tidal debris, highly asymmetric light distributions, 
double nuclei inside a common body, etc. The pres- 
ence of such distortions is considered indicative of 
a past merger that occurred between times (^obs - 
<vis) and tobs- We also make an extra test to verify 
that the distortions are caused by a past merger, 
rather than a present tidal interaction, by verifying 
that such systems do not have a distorted compan- 
ion of similar spectrophotometric redshift within 40 
kpc. 

2. They appear as a single distorted system, rather 
than 2 individually recognizable galaxies, in the 
ACS images of PSF ~0'.'l (corresponding to 380 
pc at z - 0.24 and 750 pc at z ~ 0.80). This sug- 
gests the system is an advanced merger where the 
2 progenitor galaxies have had time to coalesce into 
a single ACS system by time iobs- 

Such advanced mergers have a single redshift (Wolf et 
al. 2004), stellar mass K'h (Borch et al. 2006), and UV- 
based star formation rate (Bell et al. 2005, 2007) deter- 
mined from the COMBO-17 ground-based data of resolu- 
tion - 1'.'5. The lack of resolved COMBO-17 data for the 
individual progenitor galaxies that led to the remnant is 
not a problem because we are only concerned in the anal- 
ysis with the stellar mass and SFR of the merger. Fur- 
thermore, the condition M* > Mcut is expected to apply 
to the merger, and not to the progenitors, in both model 
and observations. 

The majority of the mergers we identify are of type Int- 
1. Examples in Figure[2]are cases 2, 3, 4, 5, 6, 7, 9, 10 11, 
and 12. 

For mergers of type Int-1 with a single redshift and stel- 
lar mass, the evidence for a merger of mass ratio > 1/10 
does not come from a measured stellar mass ratio M1/M2, 
but instead is inferred from the presence of the afore- 
mentioned morphological distortions, which are seen in 
simulations of mass ratio M1/M2 > 1/10. 



Without individual stellar masses Ml and M2, the fur- 
ther separation of these mergers into major and minor is 
not possible in every case, since the morphological dis- 
turbances induced depend not only on the mass ratio of 
the progenitors, but also on the orbital geometry (pro- 
grade or retrograde), the gas mass fraction, and structural 
parameters (e.g., Mihos & Hernquist 1996; Struck 1997; 
Naab & Burkert 2001; Mihos et al. 1995, di Matteo et 
al. 2007). Instead, we can only separate such mergers 
into three groups: clear major mergers clear minor merg- 
ers, and ambiguous cases of 'major or minor' merger as 
follows: 

1. The class of clear major mergers includes systems 
that host fairly unique tell-tale morphological dis- 
tortions characteristic of a major merger, such as 2 
tails of equal lengths, 2 nuclei of similar luminosity, 
(e.g. case 6 in Fig. [J) and a train- wreck morphology 
(e.g., case 12 in Fig. [2]). 

2. The class of clear minor merger includes the merger 
systems where the outer disk (identified based on 
the presence of disk features, such as spiral arms 
and bars) has clearly survived the recent past 
merger. Examples include case 2 of a warped disk in 
Fig. [2] The reasoning behind classifying such rem- 
nants as a minor merger is that the outer disk of 
a galaxy survives a minor merger, but not a major 
merger. While major mergers of extremely gas-rich 
disks with low star formation efficiency can lead to 
a remnant with an extended stellar disk (Robert- 
son et al. 2004), such mergers are unlikely to be 
relevant for our study, which focuses on fairly mas- 
sive systems (with stellar masses > 1 x 10^ -^o) at 
z<l. 

One further criterion is applied. When identifying 
minor mergers through the presence of a distorted 
surviving outer disk, we take care to check that the 
candidate is not in the early phases of a merger, 
which would destroy the disk in the near future, 
on a timescale ivis- We do this by avoiding those 
galaxies, which have both a distorted outer disk and 
a close companion of similar redshift (within the 
photometric redshift accuracy) and of mass ratio 
> 1/4. 

3. The class of ambiguous 'major or minor' merger is 
assigned to systems hosting distortions, which could 
be due to both a major and a minor merger. Ex- 
amples in Fig. [5] are cases 3, 4, 5, 7, 8, 10, and 
11. 

3.2.2. Mergers of type Int-2 

Mergers of type Int-2 primarily represent young merg- 
ers, which appear as very close pairs of overlapping galax- 
ies (VCPOG) in ACS images, have M* > Mcut, and are 
likely associated with a merger of mass ratio > 1/10 that 
occurred betwen times (tobs - ^vis) and iobs- 

Systems of type Int-2 are identified empirically based 
on the following criterion: ACS images show 2 individu- 
ally recognizable galaxies whose bodies overlap to form a 
common continuous envelope of light, and whose centers 
have a small separation d < 20 kpc. One or both of the 
galaxies often have morphological distortions. These prop- 
erties suggest that the 2 progenitor galaxies have recently 
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merged, at a time close to ^obs, but have not yet coalesced 
into a single ACS system of type Int-1. 

It is important to note that we are not concerned here 
with pairs of galaxies with a wide separation d > t^is x 'v 
(where v is the relative speed between the 2 progenitor 
galaxies) , because such systems represent potential future 
mergers that will occur in the next discrete time step of 
{tabs + ivis)- Instead, we are only interested in very close 
overlapping pairs of galaxies with separation d <C tvis x v, 
which represent young mergers that have occurred at a 
past time close to iobs7 and which will likely coalesce into 
a single ACS system well before the next discrete time 
step of (iobs + ivis)- For t^is of 0.5 Gyr and v ~ 100 km/s, 
d ivis X V translates to d <C 50 kpc. We thus only 
consider VCPOG where the galaxies overlap and have a 
separation d < 20 kpc, which corresponds to d < 5'.'3 at 
z - 0.24 and d < 2'.'8 at ^ ~ 0.80. 

One caveat in handling systems of type Int-2 is that 
some of the VCPOGs may be chance projections rather 
than real gravitationally bound mergers. However, this 
uncertainty does not significantly affect our result as the 
vast majority (> 80%) of the mergers in our study are of 
type Int-1 rather than Int-2. Furthermore, the likelihood 
of chance projection is mitigated due to the fact that we 
are considering pairs of very small separation. 

The COMBO-17 ground-based data of resolution - I'.'S 
will resolve a fraction of the VCPOG that make up sys- 
tems of type Int-2. Thus, we divide the latter systems into 
sub-classes Int-2a and Int-2b: 

1. Mergers of type Int-2a: These very close pairs of 
overlapping galaxies are not resolved by COMBO- 
17 data. Thus, only one stellar mass, redshift, and 
UV-based SFR are available for the pair. The lack 
of resolved COMBO-17 data for each galaxy in the 
the pair is not a problem because we are only con- 
cerned with the mass and SFR of the progenitor. 
Furthermore, the condition > Mcut is expected 
to apply to the merger, and not to the progenitors, 
in both model and observations. 

Since only one mass is available for the entire 
merged system, the evidence for a merger of mass 
ratio > 1/10 in systems of type Int-2a does not 
come from a measured mass ratio, but instead is 
deduced from the presence of morphological distor- 
tions, which are seen in simulations of mass ratio 
M1/M2 > 1/10. Essentially, the same approach 
described for systems of type Int-1 in 5 13. 2. H is used 
here. 

2. Mergers of type Int-2b: These VCPOG are resolved 
by COMBO-17 data such that stellar masses Ml 
and M2, as well as redshifts and UV-based SFRs, 
are available for both galaxies in the pair. An ex- 
ample is case 1 in Fig. [5]). For systems of type 
Int-2b, the evidence for a merger of mass ratio 
> 1/10 comes directly from the measured mass ra- 
tio M1/M2. The SFR and mass of the merger is 
considered as the sum of that of its two progenitor 
galaxies 

One caveat should be noted concerning the com- 
pleteness of mergers of type Int-2b. Strictly speak- 
ing, the criterion > Mcut applies to the merger 



mass (M1-I-M2) rather than to Ml or M2 individ- 
ually. For Mcut 2.5 x 10^" Mq counting all major 
mergers of type Int-2b with 1/4 < M1/M2 < 1/1, 
requires our sample to be complete down to 1 x 10"'^'' 
Mq for 1:1 major mergers, and down to 8 x 10^ Mq 
for 1:3 major mergers. Similarly, counting all minor 
mergers of type Int-2b with 1/10 < M1/M2 < 1/4 
requires our sample to be complete down to 6 x 10^ 
Mq for 1:4 minor mergers, and down to 2.5 x 10^ 
Mq for 1:9 minor mergers. For the blue cloud, 
where we are complete down to 1.0 x 10^ Mq, our 
count of mergers of type Int-2b is complete, but 
there will inevitably be incompleteness on the red 
sequence, particularly among minor mergers. The 
impact of this incompleteness on our results is mit- 
igated by the fact that most of the galaxies in our 
sample is on the blue cloud (see Fig. [T]) rather than 
on the red sequence, and the vast majority (> 80%) 
of the mergers in our study are of type Int-1 rather 
than Int-2. 

3.3. Visual Classification of Non-Interacting Galaxies 

Systems that do not satisfy the criteria in ij |3.2l and show 
no evidence of a recent merger of mass ratio > 1/10 are 
classified as non-interacting. These systems may harbor 
very subtle distortions, but none of the kind shown by the 
mergers of type 'Int-1 and 'Int-2'. The non-interacting 
systems are divided into two sub-groups: non-interacting 
E-to-Sd and non- interacting Irrl. These are shown in Fig- 
ure [3] and described below. 

1. 'Non-Interacting Irrl': It is important to note that 
even non-interacting galaxies have some inherent 
level of small-scale asymmetries in optical light due 
to star formation activity. In the case of low mass 
galaxies, further asymmetries may also arise due to 
the low ratio of rotational to random velocities, as 
is commonly seen in Im and Sm. These internally- 
triggered asymmetries due to SF in non-interacting 
galaxies differ in scale (few 100 pc vs several kpc) 
and morphology from the externally-triggered dis- 
tortions typical of the 'Int-1' class. We classify non- 
interacting galaxies with such internally-triggered 
asymmetries as 'Irrl' (see Figure [3l). Such systems 
may get erroneously classified as mergers in auto- 
mated asymmetry-based codes (see § 14. 2p . 

2. 'Non-Interacting E-to-Sd': Galaxies are assigned 
the 'E to Sd' class if they are fairly symmetric, have 
Hubble types in the range E-to-Sd, and are not asso- 
ciated with any overlapping or contact companion. 

In this paper, we are primarily concerned about the 
differences between three groups: the mergers in class 
'Int', the non-interacting E-to-Sd galaxies, and the non- 
interacting Irrl galaxies. The details of how E-to-Sd galax- 
ies are further sub-divided into individual Hubble types 
do not have any major impact on our main results. We 
nonetheless briefly describe this sub-classification as it is 
of interest to other studies and relevant for the test pre- 
sented at the end of § 14.11 We use conventional definitions 
(Binney & Merrifield 1998) for individual Hubble types 
(E, SO, Sa, Sb-Sc, and Sd). We assign an elliptical (E) 
type if a galaxy exhibits a smooth featureless appearance, 
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shows no disk signatures, such as a bar or spiral arms, and 
appears to be a pure spheroid. We assign an SO class if 
a galaxy hosts a smooth central brightness condensation, 
surrounded by an outer component, which is relatively fea- 
tureless (without spiral arms) and has a less steeply declin- 
ing brightness profile. We assign Sa, Sb-Sc, and Sd types 
using primarily visual estimates of the B/D ratio, and 
secondarily the smoothness/dumpiness of the disk. At in- 
termediate redshifts, where the faint smooth arms of Sa 
galaxies are not easily discernible, the distinction between 
E, SO, and Sa becomes blurred (see also ? 14. ip . However, 
this ambiguity between Es, SOs and Sas is not a problem 
for the subsequent analyses in this paper, since galaxies 
are grouped together either as 'E-|-SO-|-Sa' or 'E-to-Sd'. 

The fraction of interacting systems (i.e., mergers), 'non- 
interacting E-to-Sd' galaxies, and 'non-interacting Irrl' 
galaxies is shown in Table [1] for the high mass sample, 
and in Table [5] for the intermediate mass sample. Further 
results and tests on the merger history from visual classes 
are presented in § 14.11 and § 14.31 

3.4. CAS 

We derived the concentration C, asymmetry and 
dumpiness S (CAS) parameters by running the the CAS 
code (Conselice 2003) on the F606W images. As is stan- 
dard practice, the segmentation maps produced during the 
original source extraction (Caldwell et al. 2008) are used 
to mask neighbors on each ACS tile. The CAS code de- 
rives the asymmetry index A (Conselice 2003) by rotating 
a galaxy image by 180 deg, subtracting the rotated image 
from the original image, summing the absolute intensities 
of the residuals, and normalizing the sum to the original 
galaxy flux. CAS improves the initial input center with the 
IRAF task 'imcenter' and then performs a further refine- 
ment within a 3 X 3 grid, picking the center that minimizes 
A. The CAS concentration index C (Bershady et al. 2000) 
is proportional to the logarithm of the ratio of the 80% to 
20% curve of growth radii within 1.5 times the Petrosian 
inverted radius at r{T] — 0.2), normalized by a logarithm 

C = 5 X Iog(r8o%/r2o%) (2) 

The dumpiness index S (Conselice 2003) is defined as 
the ratio of the amount of light contained in high-frequency 
structures to the total amount of light in the galaxy. In or- 
der to compute S, the CAS code first smooths the galaxy 
image with a filter of size equal to 1/6 of the Petrosian 
radius to produce a lower resolution image whose high- 
frequency structure has been washed out. The latter im- 
age is then subtracted from the original image to produce 
a residual map that contains only the high-frequency com- 
ponents of the galaxy's stellar light. The flux of this resid- 
ual light is then summed and divided by the sum of the 
original galaxy image flux to obtain a galaxy's dumpiness 
{S) value. Tests on the interaction history from CAS are 
presented in § 14.21 

It has been argued that the criterion A > 0.35 and 
A > S (henceforth referred to as the CAS merger criterion) 
captures galaxies that exhibit large asymmetries produced 
by major mergers (Conselice 2003). We will assess this in 
§1121 

4. RESULTS AND DISCUSSION 



4.1. The merger fraction from visual classes 

Fig m compares the fraction / of systems with evidence 
of a recent merger of mass ratio > 1/10, based on vi- 
sual classification by 3 classifiers (SJ, SM, KP). Results 
are shown for both the high mass (M > 2.5 x 10^° Mq) 
sample SI and the intermediate mass (M > 1 x 10^ Mq) 
sample S2. On this figure, the plotted error bar for the 
merger fraction / only includes the binomial term [/(l- 
f)/N]^^^, for each bin of size N. The same trend is seen 
for all 3 classifiers and the maximum spread Sf/f in the 
four bins is 15%, 17%, 26%, and 26%, respectively. In 
subsequent analyses, we adopt a conservative error bar on 
/ that includes in quadrature both the binomial term and 
a dispersion of 26% to capture the inherent subjectivity in 
the visual classification. 

Another key test is to assess the impact of redshift- 
dependent systematic effects, such as bandpass shifting. 
When using the F606W filter whose pivot wavelength is 
~ 5915 A, the rest frame wavelength (Arcst) corresponds 
to the rest-frame optical at the mean redshift of the first 
3 bins, but shifts to the rest-frame violet /near- UV (3700 
A to 3290 A) in the last bin (z ^ 0.6 to 0.8). Galaxies 
tend to look slightly more asymmetric at near-UV wave- 
lengths due to the prominence of young stars. In order 
to quantitatively test the impact of bandpass shift on our 
visual classes, we use the redder F850LP images from the 
GOODS survey, which overlaps with the central 20% of 
the GEMS survey area. The F850LP filter has a pivot 
wavelength of 9103 A and traces the rest-frame optical 
(7340 A to 5057 A) in all four redshift bins out to z ~ 0.8. 
The F850LP images also have 5 times longer exposures 
than the GEMS F850LP and F606W images. Figure [5] 
shows GEMS F606W and GOODS F850LP images of typ- 
ical disturbed and normal galaxies in the last 2 redshift 
bins (z ~ 0.47 to 0.8). 

While the GOODS images have higher S/N, and trace 
redder older stars, they do not reveal dramatically differ- 
ent morphologies from those in the GEMS F606W im- 
ages (Fig [SJ. Furthermore, the 855 intermediate mass 
{M > 1 X 10^ Mq) galaxies in the GEMS/GOODS over- 
lap area, were classified using both GOODS F850LP and 
GEMS F606W images by the 3 classifiers. We find that the 
ratio of (/gems//goods) ranges from 0.8 to 1.2 across the 
3 classifiers (Table [3]), where /gems and /goods are the 
fraction of merging systems based on the GEMS F606W 
and GOODS F850LP images, respectively. The mean / 
changes by only 6% (Table [3]). In effect, over 85% of the 
systems classified as mergers ('Int') in the GEMS F606W 
images retain the same visual class in the GOODS F850LP. 
Among the remaining objects, some classified as non- 
interacting in GEMS F606W get reclassified as disturbed 
in GOODS F850LP, and vice-versa. The fact that / does 
not change by a large amount between GEMS F606W and 
GOODS F850LP is not surprising, since the rest-frame 
wavelength of GEMS F606W in the last bin shifts only to 
the violet /near-UV, rather than to the far-UV, where mor- 
phological changes are more dramatic. We conclude that 
our results are not highly impacted by bandpass shifting, 
and any effect is accounted for by our error bars of > 26% 
in /. 

Another redshift-dependent systematic effect is surface 
brightness dimming at the rate of (l-|-z)~^ (e.g., Harden 
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et al. 2008). This leads to surface brightness dimming by a 
factor of 1.0 to 2.5 magnitude over the redshift range 0.24 
to 0.80. This is mitigated in part by two factors: galaxies 
are on average 1.0 magnitude brighter in surface brightness 
by z ~ 0.8 (e.g., Barden et al. 2005), and the average SFR 
rises by a factor of ~ 4 out to z 0.8 (e.g., see ij |4.6p . Two 
approaches can be adopted to assess the impact of surface 
brightness dimming. The first is to artificially redshift dis- 
turbed galaxies in the lowest redshift bin (z ^ 0.24) out 
to z 0.8, either assuming passive evolution or adding 
in a ~ 1 magnitude of brightening in surface brightness. 
However, this approach suffers from the limitation that it 
implicitly assumes that galaxies at z ~ 0.8 are similar to 
those at z ~ 0.24 and evolve passively with time. A better 
approach, which does not make such assumptions, is to 
repeat the analysis and visual classification using deeper 
images of the galaxies and assess the resulting change in 
visual classes. The above-described test performed using 
the deep GOODS F850LP image (Fig. is an example 
of such a test, and indicates that the eye-ball morpholo- 
gies do not change within the error bars of > 26% in /. 
We note however that quantitative CAS parameters can 
change with the deeper GOODS images (e.g., Conselice et 
al. 2008). 

Finally, as an extra test, we checked the distribution 
of Sersic indices n for single-component Sersic fits (Bar- 
den et al. 2005) for the visual classes of the sample S2 of 
intermediate mass M > 1 x 10^ Mq galaxies (Fig. (5]). Non- 
interacting disk-dominated systems are expected to have 
n < 2.5, while massive ellipticals and bulge-dominated 
systems typically have higher Sersic indices. We indeed 
find that over 85% of the systems visually classified as 
Sb-Sd and Irrl have n < 2.5 in the intermediate mass 
(M > 1 X 10^ -^o) sample. Furthermore, as expected, 
the vast majority of the systems typed as Sa have n < 4. 
However, the systems typed as E and SO span a broad 
range in n: most of them have n > 3, but there is a tail 
of lower n values. This is not surprising given the previ- 
ously described difficulties fij l3.3p in separating E, SO, and 
Sa galaxies at intermediate redshifts. However, this am- 
biguity between E, SO, and Sa systems is not a problem 
for the subsequent analyses in this paper, since galaxies 
are grouped together either as 'E-|-SO-|-Sa' or 'E-to-Sd'. In 
fact, as stressed in § 13. 3[ the main results presented in 
this paper depend only on the differences between three 
groups: mergers (Tnf), non-Interacting E to Sd galaxies, 
and non-interacting Irrl galaxies. 

4.2. The merger fraction from CAS 

It has been argued that the CAS merger criterion 
{A > 0.35 and A > S) captures systems that exhibit 
large asymmetries produced by major mergers (Conselice 
2003). This criterion is based on calibrations of the CAS 
system at optical rest-frame wavelengths (Arest > 4500 A) . 
However, there are several caveats: a) The CAS criterion 
{A > 0.35 and A > S) will miss out interacting galaxies 
where the morphological distortions contribute to less than 
35 % of the total galaxy flux, (b) Calibrations of A with 
N-body simulations (Conselice 2006) shows that during 
major mergers with mass ratios 1:1 to 1:3, the asymmetry 
oscillates with time. Typically, it exceeds 0.35 for ^ 0.2 
Gyr in the early phases when the galaxies start to interact. 



falls to low values as the galaxies separate, rises for ~ 0.2 
Gyr as they approach again for the final merger, and even- 
tually tapers down as the final remnant relaxes. On aver- 
age, the A > 0.35 criterion is only satisfied for one third of 
the merger timescale in these N-body simulations. For mi- 
nor mergers of mass ratios 1:5 and below, the asymmetries 
are too low to satisfy A > 0.35. (c) To complicate matters, 
automated asymmetry parameters can also capture non- 
interacting galaxies whose visible light shows small-scale 
asymmetries due to star formation (e.g.. Miller et al. 2008; 
Lotz et al. 2008). 

Visual tests that verify how well the CAS criterion 
(A > 0.35 and A > S) works at intermediate redshifts have 
been performed using spot checks and small-to-moderate 
samples (e.g., Mobasher et al 2004; Conselice 2003; Con- 
selice et al. 2003; Conselice et al. 2005). However, what 
has been missing to date is a quantitative estimate, based 
on a large sample of galaxies, of the recovery fraction of 
CAS (i.e., the fraction of visually-classified mergers that 
the CAS criterion picks up) , and the contamination level of 
CAS (i.e., the fraction of visually-classified non-interacting 
galaxies that the CAS criterion picks up). Both the recov- 
ery fraction and contamination level might be expected to 
depend on the rest-frame wavelength used, the mass and 
SFR of the galaxies, etc. In this paper, we perform one of 
the most extensive comparisons to date, at intermediate 
redshifts (z ~ 0.24 to 0.80), between CAS-based and visual 
classification results for both high mass (M > 2.5 x 10^" 
Mq) and intermediate mass (M > 1 x 10^ Mq) galaxies. 
We assess the effectiveness of the CAS merger criterion 
{A > 0.35 and A > S) over this interval, where the rest- 
frame wavelength Arost varies from 4770 A to 3286 A. We 
note that the rest-frame wavelength range here extends to 
somewhat bluer wavelengths than the range (Arost > 4500 
A) over which the CAS system was calibrated. 

Fig |4] compares the merger fractions that would be ob- 
tained using the CAS criterion (/cas), as opposed to visual 
classification (/). For the high mass {M > 2.5 x 10^° Mq) 
galaxies, visually based and CAS-based merger fractions 
agree within a factor of two, with / being higher than /cas 
at z < 0.5, and being lower at z > 0.5 (top panel of Fig|4]). 
However, for the intermediate mass (M > 1 x 10^ Mq) 
galaxies (lower panel of FigH]), at z > 0.5 the CAS-based 
merger fraction can be systematically higher by a factor 
~ 3 than the visually based /. The reason for this discrep- 
ancy, as we show below, is that at bluer rest-frame wave- 
lengths (i.e., higher redshifts), the CAS criterion picks up 
a significant number of non-interacting dusty, star-forming 
galaxies. 

Fig. [7| plots the CAS asymmetry A and dumpiness S 
parameter for galaxies in the four redshift bins covering 
the interval z ~ 0.24-0.80. Galaxies satisfying the CAS 
criterion {A > 0.35 and A > S) lie in the upper left hand 
corner. One can see that while the CAS criterion captures 
a fair fraction of the mergers (coded as orange stars), it 
also picks up a large number of non-interacting galaxies. 

We define the recovery fraction (FcAS-merger) of CAS 
as the fraction of visually-classified mergers ('Int'), which 
are picked up by the CAS criterion {A > 0.35 and A > S). 
For the high mass sample (M > 2.5 x 10^" Mq) sample, 
the recovery fractions in the four redshift bins are 50% 
(2/4), 14% (1/7), 42% (8/19), and 56% (20/36) respec- 
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tively, with low number statistics dominating the first two 
bins. For the intermediate mass (M > 1 x 10^ -^0) sam- 
ple, the recovery fractions in the four redshift bins are 50% 
(13/26), 69% (30/43), 59% (49/83), and 73% (85//116), 
respectively, as illustrated in the top panel of Fig. \8\ We 
inspected the visually-classified mergers missed out by the 
CAS criterion {A > 0.35 and A > S) and show typical 
cases in the top panel of Fig. [9l The missed cases include 
galaxies where tidal or accretion features in the main disk 
of a galaxy contribute less than 35% of the total hght (e.g., 
case 3 in Fig. [HI); galaxies with close double nuclei (e.g., 
case 2 in Fig. [5]) where CAS might refine the center to be 
between the two nuclei, thereby leading to a low A < 0.35; 
and pairs of fairly symmetric galaxies whose members have 
similar redshifts within the spectrophotometric error, ap- 
pear connected via weak tidal features, and have a stellar 
mass ratio Afl/M2 > 1/10 (e.g., case 1 in Fig. [5] where 
M1/M2 ~ 0.25). 

We define the contamination fraction of CAS as the 
fraction of those systems which satisfy the CAS crite- 
rion {A > 0.35 and A > S) and are therefore consid- 
ered as likely major mergers by CAS, but yet are visually 
classified as non-interacting. For the high mass sample 
(M > 2.5 X 10^° Mq) sample, the CAS contamination 
fractions in the four redshift bins are 34%, 75%, 72%, and 
67% respectively, with low number statistics dominating 
the first two bins. For the intermediate mass (M > 1 x 10^ 
Mq) sample, the corresponding CAS contamination frac- 
tions are 44%, 53%, 76%, and 82% respectively, as shown 
in the lower panel of Fig. [51 On the latter figure, Ncas 
represents the total number of galaxies satisfying the CAS 
criterion {A > 0.35 and A > S) in the four redshift bins. 
Plotted on the y-axis is the fraction i^cAS-visuai of differ- 
ent visual types (mergers, non-interacting E-Sd, and non- 
interacting Irrl) among these "CAS mergers". Across the 
four redshift bins, the non-interacting E-Sd and Irrl make 
up 44%, 53%, 76%, and 82% of the CAS systems. Typical 
cases are shown in the lower panel of Fig. [51 They in- 
clude non-interacting actively star-forming systems where 
SF induces small-scale asymmetries in the optical blue 
light (e.g., cases 4 and 6 in Fig. [5]); systems where A is 
high due to the absence of a clearly defined center (e.g., 
case 8 in Fig. [9]) or due to the center being blocked by 
dust (e.g., cases 4 and 9 in Fig. [S]); and compact or edge- 
on systems where the light profile is steep such that small 
centering inaccuracies can lead to large A (e.g., case 9 in 
Fig. [ID. 

In summary, we find that the CAS-based merger frac- 
tion agrees within a factor of two with visually based one 
for high mass (M > 2.5 x 10^" Mq) galaxies, but can over- 
estimate the merger fraction at z > 0.5 by a factor 3 for 
intermediate mass (M > 1 x 10^ Mq) galaxies. For the 
latter mass range, the systems counting toward /cas are 
a mixed bag: the CAS criterion misses about half of the 
visually-classified mergers, but picks up a dominant num- 
ber of non- interacting dusty, star- forming galaxies. We 
thus conclude that the CAS merger criterion is ill-suited 
for use on HST F-band images at z > 0.5, where the rest 
frame wavelength falls below A < 4000 A, particularly in 
the case of intermediate mass galaxies with significant SF, 
gas, and dust. Modified CAS criteria in the near-UV based 
on morphological k-corrections (Taylor et al. 2007) might 



alleviate this problem. 

4.3. Interaction history of massive and intermediate 
mass galaxies 

Based on the tests in 14. H and W.'2l we decided to adopt 
the mean merger fraction / based on visual classes for our 
two samples of interest. For the high mass (M > 2.5 x 10^° 
Mq) sample, which is complete on both the blue cloud and 
red sequence (§[2), the results are shown on in Table[Tl The 
error bar shown on the merger fraction / in both tables 
now includes the sum in quadrature of a binomial term 
[f{l-f)/N]^^^ for each bin of size N, along with a frac- 
tional error of ± 26% to capture the dispersion between 
classifiers, and uncertainties due to bandpass shifting and 
surface brightness dimming. 

From Tableland FiglH it can be seen that the merger 
fraction / among high mass (M > 2.5 x 10^° -^o) galax- 
ies does not show strong evolution over lookback times of 
3-7 Gyr, ranging from 9% ± 5% at z ~ 0.24-0.34, to 8% 
± 2% at z ^ 0.60-0.80, as averaged over every Gyr bin. As 
discussed in § 13. 2[ the merger fraction / refers to systems 
with evidence of a recent merger of mass ratio > 1/10. 

As outlined in § 13.21 these mergers were further subdi- 
vided among 3 classes: clear major merger, clear minor 
merger, and ambiguous 'major or minor merger' cases. 
The first two classes are used to set lower limits on the 
major and minor merger fraction. The lower limit on the 
major (Ml/Af2 > 1/4) merger fraction, determined in this 
way, ranges from 1.1% to 3.5% over z - 0.24-0.80. The 
corresponding lower limit on the minor (1/10 < Ml /M 2 < 
1/4) merger fraction ranges from 3.6% to 7.5%. To our 
knowledge, this is the first, albeit approximate, empirical 
estimate of the frequency of minor mergers over the last 
7 Gyr. The ambiguous cases of 'major or minor merger' 
make up a fraction between 1.2% to 2.0%. 

When converting the observed fraction / of galaxy merg- 
ers into a merger rate R, we must bear in mind that in 
any observational survey of galaxies, mergers can only be 
recognized for a finite time tvis, which is the timescale 
over which a merging galaxy will appear morphologically 
distorted. This timescale depends on the mass ratio of 
the merger as well as the gas fraction of the progenitors: 
ivis 0.5-0.8 for gas-rich galaxies, and tvis ~ 0.2-0.4 Gyr 
for gas-poor galaxies (T.J. Cox, private communication). 
This timescale will also depend on many observational fac- 
tors such as the method used to identify mergers (e.g. vi- 
sual classification vs. CAS or other statistical methods) 
and the depth of the imaging used. We assume a repre- 
sentative value of ivis = 0.5 Gyr here, but we must keep 
in mind that there are at least factors of two uncertainty 
in this number. The merger rate R is given by 

R=}^, (3) 

^vis 

where n is the comoving number density of galaxies above 
a certain mass limit in the redshift bin. 

For the sample of high mass (M > 2.5 x 10^° Mq) galax- 
ies, our measured merger fraction / and assumed value of 
tvis ^ 0.5 Gyr lead to a corresponding merger rate i? of a 
few xlO"'' galaxies Gyr~^ Mpc~'^. Assuming a visibility 
timescale of ^ 0.5 Gyr, it follows that on average, over 
Th 5-7 Gyr, ^ 68% of high mass systems have under- 
gone a merger of mass ratio > 1/10. Of these, we estimate 
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that ~ 16%, 4-5%, and 7% correspond respectively to clear 
major mergers, clear minor mergers, and ambiguous cases 
of 'major or minor' mergers. 

At intermediate masses (M > 1 x 10^ ^-^0) where we 
are only complete in mass for the blue cloud (§[11), we con- 
sider / to be meaningful only for the intermediate mass 
blue cloud sample. Results for this sample are shown in 
the lower part of Tabled The fraction of blue cloud galax- 
ies having undergone recent mergers of mass ratio > 1/10 
ranges from 7% ± 2% to 15% ± 5% over z - 0.24-0.80. 
The corresponding merger rate R ranges from 8 x lO"** to 
1 X 10"'^ galaxies Gyr^^ Mpc^^ . For an assumed visibil- 
ity time of ^ 0.5 Gyr, we estimate that on average, over 
Tb ~ 3-7 Gyr, 84% of intermediate mass blue cloud galax- 
ies have undergone a merger of mass ratio > 1/10, with 

5%, 22%, and 57% corresponding respectively to clear 
major mergers, clear minor mergers, and ambiguous cases 
of 'major or minor' mergers 

4.4. Comparison with other studies 

When comparing our observed merger fraction / in the 
high mass [M > 2.5 x 10^" Mq) sample over z ~ 0.24-0.80 
with published studies, several caveats must be borne in 
mind. Many studies have small samples and large error 
bars at z < 0.8 (e.g., Consehce 2003; Fig. [TUl). Others 
focus on bright galaxies and luminosity-selected samples 
(e.g., Lotz et al. 2008; Casatta et al. 2005 ) rather than 
stellar mass selected sample, because the data to derive 
stellar masses were unavailable. Different studies target 
different systems, ranging from morphologically distorted 
systems to close pairs with separation d ~ 5 to 40 kpc. 
Finally, many studies focus only on major mergers, while 
the interacting galaxies identified in our study are candi- 
dates for a merger of mass ratio > 1/10 (i) l3.2p . and include 
both minor and major mergers. Nonetheless, we attempt 
approximate comparisons. 

Fig. [TO] shows the merger fraction based primarily on 
morphologically distorted galaxies (filled circles), as well 
as the close pair fraction (open squares), as a function 
of redshift. The Lotz et al. (2008) study shows the frac- 
tion /cini of morphologically disturbed systems based on 
Gini-M20 parameters among Mb < -20.5 and Lb > 0.4 
galaxies in the Extended Groth Strip. This study does not 
present any results for a high mass sample, and thus we 
effectively are comparing their bright galaxies to our high 
mass galaxies. Over z ~0.2-0.80, our results are in very 
good agreement, within a factor of less than two, with 
/cini- The CAS-based results from Conselice (2003) are 
derived from a small sample in the Hubble Deep Field and 
have error bars that are too large to set useful constraints 
at 2; < 1 (Fig. [TO]) . Our results of a fairly flat evolution of 
the merger rate out to z ^ 0.8 also agree with the results 
of Cassata et al. (2005), which are based on both pairs 
and asymmetries. 

Fig. [TO] also shows the result from three studies based on 
close pairs. The major merger fraction of massive galax- 
ies (M, > 2.5 X 10^° Mq) in close {d < 30 kpc) pairs, 
based on the 2-point correlation function in COMBO-17, 
is 5% ± 1% averaged over at 0.4 < z < 0.8 (Bell et 
al. 2006). This value is lower than our merger fraction 
/ (~ 8% ± 2%), which represent likely mergers of mass 
ratio > 1/10, and it is higher than the fraction of cases we 



see as clear major mergers 1.3% ± 0.2%). The study 
of luminous {Ly > 0.4 L«) pairs at projected separations 
of 5-20 kpc in the COSMOS field (Kartaltepe et al. 2007) 
finds a galaxy pair fraction of ~ l%-3% over z ~ 0.24- 
0.80, corresponding to a galaxy merger fraction of ~ 2%- 
6% . Our observed fraction / of 9% ± 5 to 8% ± 2% over 
z ~ 0.24-0.8 is slightly higher and flatter than this study. 

The differences we see through these comparisons are 
already known (see § [1]) . Studies based on close pairs tend 
to show moderate to fairly strong evolution in the major 
merger rate out to z '-^ 1.2 (e.g., Kartaltepe et al 2007; Bell 
et al 2006; Lin 2004), while studies based on asymmetries 
(e.g., Lotz et al. 2008; this study), and studies based on 
both pairs and asymmetries (Cassata et al. 2005) tend to 
report only mild evolution of the merger rate with redshift 
up to z~ 1. 

It is not fully understood why these different methods 
yield different results, but several factors likely play a part. 
First, it should be noted that the claim of strong evolution 
in the close pair fraction out to z ^ 1.2 in the COSMOS 
study by Kartaltepe et al. (2007) comes about when the 
low redshift z ~ point from the SDSS pair catalog (Al- 
1am et al. 2004) is included in their analysis. The evo- 
lution within the internally consistent dataset from COS- 
MOS over z ~ 0.15 to 1.05 shows much weaker evolution 
(Fig. I10(). The drop in close pair fraction seems primarily 
to occur at z < 0.2, but it is unclear how reliable the low 
z ~ points are due to the small volume sampled and sys- 
tematic effects between studies. A further reason for the 
difference could be due to the fact that the methods used in 
these studies trace different phases of an interaction, with 
the pair method tracing the potential pre-merger phase, 
while the method based on the distorted galaxies trace 
the later phases, including the merger and post-merger 
phases. 

Another point is that both pair and asymmetry meth- 
ods are imperfect ways of tracing the merger fraction. 
Methods tracing morphologically disturbed galaxies may 
capture some fly-by tidal interactions rather than merg- 
ers, and this effect would cause the fraction of interacting 
galaxies to overestimate the merger fraction. However, 
this effect is not a dominant one due to the following rea- 
son: interaction signatures typically persist for a visibility 
timescale of 0.5 Gyr (Tvis), and a fly- by companion caus- 
ing the distortion would still be within 100 kpc of the 
disturbed galaxy, assuming an escape speed of 200 km/s. 
The distorted galaxies we identify do not typically have 
such a fly- by companion, of mass ratio > 1/10 and simi- 
lar spectrophotometric redshift. In studies based on close 
pairs, one source of uncertainty is that even pairs with 
members at the same redshift may not become gravitation- 
ally bound in the future. This effect might cause pairs to 
overestimate the true major merger fraction. On the other 
hand, erroneous spectrophotometric redshifts can cause us 
to either overestimate or underestimate the true close pair 
fraction, with the latter effect being more likely. Correc- 
tions for this effect are uncertain and depend on the shape 
of the spectrophotometric redshift errors (e.g., see Bell et 
al. 2006 for discussion). 

4.5. Comparison of galaxy merger history with ACDM 

models 
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We compare our empirical merger fraction / (Fig. Ilip 
and merger rate R (Fig. [T2)) to predictions from different 
theoretical models of galaxy evolution in the context of a 
ACDM cosmology, including the halo occupation distribu- 
tion (HOD) models of Hopkins et al. (2007); semi-analytic 
models (SAMs) of Somerville et al. (2008), Bower et al. 
(2006), and Khochfar & Silk (2006); and the cosmologi- 
cal smoothed particle hydrodynamics (SPH) simulations 
from Mailer et al. (2006). The models were provided to us 
directly by the authors or co-authors of these individual 
studies. 

We first briefly describe the general problem of calcu- 
lating galaxy merger rates. Predicting the rate of merg- 
ers per comoving volume and per unit time between iso- 
lated dark matter (DM) halos within a ACDM model is 
relatively straightforward via semi-analytic methods or N- 
body simulations (e.g. Lacey & Cole 1993; Gottlober et 
al. 2001; Fakhouri & Ma 2008; Neistein & Dekel 2008; 
D'Onghia et al. 2008). However, making a direct predic- 
tion of the galaxy merger rates is more complicated due 
to a number of factors, including the difference between 
the galaxy and halo merger timescales, tidal heating and 
stripping of halos and sub-halos, the effect of a dense core 
of baryons on merging satellites, and the non-linear rela- 
tion at low mass between DM halo (or sub-halo) mass and 
galaxy mass (van den Bosch et al. 2007). Thus, attempts 
to extract a galaxy merger rate from ACDM simulations 
also must attempt to model the relationship between dark 
matter and galaxy properties. The three main methods for 
making this connection are HOD models, SAMs, and hy- 
drodynamic simulations. We summarize below how these 
three types of models differ. 

HOD models specify the probability that a DM halo of 
a given mass M harbors N galaxies above a given mass 
or luminosity. The parameters of this function are de- 
termined by requiring that statistical observed quantities, 
such as galaxy mass or luminosity functions and galaxy 
correlation functions, be reproduced. The merger rate of 
galaxies within their host halos is calculated via standard 
or improved dynamical friction formulae. In the HOD 
models of Hopkins et al. (2007) used here, different mod- 
ified formulae can be used, which include the effect of a 
gravitational capture cross section, stripping of DM ha- 
los, and calibration factors from A'^-body simulations. The 
predicted model rate can vary by a factor of ~ two de- 
pending on model assumptions for sub-halo structure and 
mass functions, the halo occupation statistics, and the dy- 
namical friction formulae used. 

In SAMS, merger trees of DM halos are either extracted 
from cosmological N-body simulations or derived using 
analytic methods (e.g., Somerville & Kolatt 1999). Cal- 
ibrated modified versions of the Chandrasekhar dynam- 
ical friction approximation (e.g. Boylan-Kolchin et al. 
2008) are used to compute the galaxy merger rate. Sim- 
plified analytic formulae are used to model the cooling 
of gas, star formation, supernova feedback, and more re- 
cently, AGN feedback (e.g., Somerville et al. 2008; Bower 
et al. 2006; Croton et al. 2006; Benson et al. 2005; Cole 
et al. 2000; Somerville & Primack 1999). The free pa- 
rameters in these formulae are normalized to reproduce 
observations of nearby galaxies, such as the z = galaxy 



mass or luminosity function. Fig. ll2l shows results of three 
independent SAMs from Khochfar & Silk (2006), Bower et 
al. (2006), and Somerville et al. (2008). 

Cosmological hydrodynamic simulations attempt to 
model the detailed physics of gas hydrodynamics and cool- 
ing as well as gravity by explicitly solving the relevant 
equations for particles or grid cells. SPH methods are 
most commonly used. SF and supernova feedback are 
treated using empirical recipes. A drawback of this ap- 
proach is that, due to computational limitations, state-of- 
the-art simulations still do not have the dynamic range to 
resolve the internal structure of galaxies while simultane- 
ously treating representative cosmological volumes. An- 
other well known problem is that cosmological SPH mod- 
els, which do not include some kind of suppression of cool- 
ing (e.g., due to AGN feedback) in massive halos do not 
reproduce the observed number density of galaxies on the 
mass scales of interest (few x10^°Mq). Thus, the simula- 
tions of Mailer et al. (2006) shown here, over-predict the 
number of high and low mass galaxies, while galaxies at 
the bend of the Schechter mass function are a factor of 2 to 
3 too massive. In order to make the simulated mass func- 
tion agree better with observations. Mailer et al. (2006) 
apply a correction factor of 2.75 for galaxies in the mass 
range 2 x 10^" < M^/Mq < 6 x 10". This correction is 
already included in the model on Fig. [T2l 

When comparing the observations to the models, one 
must consider carefully how merger rates and fractions 
are determined in these simulations. Two approaches are 
used: one based on simulation snapshots and the other 
based on a light cone. In the first approach, simulation 
outputs ( "snapshots" ) are stored at a sequence of redshifts. 
Two snapshots separated by a time At are considered and 
modelers trace the merger history of galaxies whose final 
stellar mass is greater or equal to a given mass cut 
Mcut- The same mass cut (M* > 2.5 x 10^°) used in the 
data is applied to the simulations. In order to mimick the 
observations as closely as possible, the interval At in the 
model should ideally be equal to the visibility timescale 
^vii3- One then counts the number iVi of model galaxies 
with M^, > Mcut, which have experienced a merger of mass 
ratio M1/M2 > 1/10 within the last tvis- It is important 
to note two points. Firstly, if a galaxy were to undergo 
multiple mergers within a time ivis, these mergers would be 
counted only once in the model term A^i , analogous to the 
case of the data where one cannot discriminate between 
multiple mergers within the time tvis over which a galaxy 
appears distorted due to a merger. Secondly, the fact that 
tvis is so short makes it very unlikely for a model galaxy 
to undergo more than one merger over this timescale. As 
a result the model term A^i is essentially tracing the num- 
ber A'mrgi of mergers within the last tvis- One can then 
determine the merger rate using i?modi = A^mrgi/(At V), 
where V is the comoving volume of the simulation box. 
The merger fraction is /modi = A'^mrgi/A'^gaii, where iVgaii 
is the total number of galaxies above the relevant mass 
limit. Except for Somerville et al. (2008), all the models 
presented on Fig. [TT] and FiglT^] derive the merger fraction 
/ and rate R using the above approach, based on simula- 
tion snapshots. 

For the Somerville et al. (2008) models, the simulation 



If At is larger than tyis, a correction factor of order (tvis/At) needs to be applied to the model merger fraction /modi 
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analysis was carried out in a way that is closer to the 
observations. We construct a light cone with a geome- 
try that is equivalent to three GEMS fields (2700 arcmin^ 
from 0.1 < z < 1.1). We then divide the galaxies into 
redshift bins, exactly as in the observational analysis, and 
we count the number N^,-g2 of galaxies that have had a 
merger within a time tvis in the past. These galaxies ap- 
pear as morphologically distorted mergers in the observa- 
tions. The model merger fraction is /mod2 = -^mrg2 /^gai2 7 
where A^gai2 is the total number of galaxies above the rel- 
evant mass limit in this light cone or redshfit bin. The 
model merger rate is calculated exactly as in the data us- 
ing i?niod2 = /mod2 rimod2/tvis, wherc nmod2 is the comov- 
ing number density of galaxies above a certain mass limit 
in the redshift bin. Note that A^mrg2 and /mod2 are quite 
sensitive to tvis, while R is independent of ivis- 

In Fig. 1111 we compare the empirical merger fraction / 
to the corresponding model predictions /modi and /mod2- 
The comparison between the empirical merger rate R 
(equation (3) in § 14. 3p and the model predictions i?modi 
and i?mod2 is in Fig. 1121 In both data and models, major 
and minor mergers are defined as those with mass ratio 
(M1/M2 > 1/4), and (1/10 < M1/M2 < 1/4), respec- 
tively. The only slight exception is in the case of Mailer et 
al (2006) model where the extracted major mergers were 
defined with a slightly lower mass cutoff (Ml/Af2 > 1/3). 
The dotted lines on Fig. [TT] and Fig. [T^ show the major 
merger rate for all the models. The solid line show the 
( maj or -|- minor) merger rate, in other words, the rate of 
mergers with mass ratio (M1/M2 > 1/10). This is shown 
for all models except the Mailer et al (2006) SPH sim- 
ulations, where the limited dynamic range of the current 
simulations only allows predictions for major mergers. The 
solid line can be directly compared to our empirical merger 
fraction / or rate R, since the latter refer to mergers of 
mass ratio M1/M2 > 1/10. 

We find qualitative agreement between the observations 
and models, such that the (major+minor) merger fraction 
(Fig. and merger rate (Fig. from different mod- 
els (solid lines) bracket the corresponding empirical esti- 
mates (stars) and show a factor of five dispersion. One 
can now anticipate that in the near future, improvements 
in both the observational estimates and model predictions 
will start to rule out certain merger scenarios and refine 
our understanding of the merger history of galaxies. 

4.6. The impact of galaxy mergers on the average SFR 
over the last 7 Gyr 

Both observations (e.g., Larson & Tinsley 1978; Joseph 
& Wright 1985; Kennicutt et al. 1987; Barton et al 2003) 
and simulations (e.g., Negroponte & White 1983; Hern- 
quist 1989; Barnes & Hernquist 1991, 1996; Mihos & Hern- 
quist 1994, 1996; Springel, Di Matteo & Hernquist 2005b) 
suggest that galaxy interactions and mergers trigger star 
formation. However, simulations cannot uniquely predict 
the factor by which mergers enhance the SF activity of 
galaxies over the last 7 Gyr, since the star formation rate 
in simulations is highly sensitive to the stellar feedback 
model, the bulge-to-disk (B/D) ratio, the gas mass frac- 
tions, and orbital geometry (e.g.. Cox et al 2006; di Matteo 
et al. 2007). Thus, we explore here the impact of inter- 
actions on the average UV-based and UV-|-IR-based SFR 



of intermediate-to-high mass (A/ > 1 x 10^ Mq) galaxies 
over z - 0.24-0.80. 

We adopt the SFRs in BeU et al. (2005, 2007), based on 
COMBO-17 UV data (Wolf et al. 2004) and deep Spitzer 
24 /im observations with a limiting flux of ^ 83 Jy (5cr) 
from the Spitzer Guaranteed Time Observers (Papovich et 
al. 2004; Gordon et al. 2005). The unobscured SFR based 
on the directly observable UV light from young stars was 
computed using SFRuv = 9.8 x 10^^^ (2.2 £uv), where 
Luv = 1-5j^L,2800 is a rough estimate of the total inte- 
grated 1216-3000 A UV luminosity, derived using the 2800 
A rest-frame luminosity from COMBO-17 ^!y,2800- The fac- 
tor of 1.5 used in converting the 2800 A luminosity to to- 
tal UV luminosity accounts for the UV spectral shape of 
a 100 Myr-old population with constant SFR. The factor 
of 2.2 corrects for the light emitted longward of 3000 A 
and shortward of 1216 A. The SFR calibration is derived 
from Pegase assuming a 100 Myr old stellar population 
with constant SFR and a Chabrier (2003) IMF. 

The obscured SFR can be calculated from 
dust-reprocessed IR emission using the expression 
SFRiR = 9.8 X 10"" XiR, where Lir is the total IR 
luminosity (TIR) over 8-1000 ^m (Bell et al. 2007). 
Lir is constructed from the observed 24 /im flux (cor- 
responding to rest-frame wavelengths of 19-13 fiui over 
z ^ 0.24-0.80) using the method outhned in Papovich & 
Bell (2002), based on an average Sbc template from the 
Devriendt et al. (1999) SED library. In converting from 
Lir to SFRiR, Bell et al. (2007) assume that the bulk of 
the 24 /im emission comes from SF, and not from AGN 
activity, based on the statistical result that less than 15% 
of the total 24 /im emission at z < 1 is in X-ray luminous 
AGN (e.g., Silva et al. 2004; BeU et al. 2005; Franccschini 
et al. 2005; Brand et al. 2006). Uncertainties in these 
SFR estimates are no less than a factor of 2 for individual 
galaxies while the systematic uncertainty in the overall 
SFR scale is likely to be less than a factor of 2 (Bell et 
al. 2007). 

We investigate the star formation properties of the sam- 
ple SI of - 789 high mass (M > 2.5 x lO^" Mq) and the 
sample S2 of ^--^3698 intermediate mass {M > 1.0 x 10^ 
Mq) galaxies. As described in §[2l the high mass sample 
SI is complete for both the red sequence and blue cloud, 
while the intermediate mass sample S2 is only complete 
for the blue cloud and suffers from incompleteness on the 
red sequence in the highest redshift bins. However, since 
most of the SFR density originates from the blue cloud, 
this incompleteness does not have any major impact on the 
results. Fig. [T3] shows the UV-based SFR plotted versus 
the stellar mass in each redshift bin. The UV-based SFR 
ranges from 0.01 to 25 Mq yr~^, with most galaxies 
having a rate below 5 Mq yr~^. 

While it is desirable to use the Spitzer 24 /i,m data in 
order to account for obscured star formation, only ^ 24% 
878 galaxies) of the 3698 galaxies in our intermediate 
mass sample have a Spitzer 24 /im detection, although over 
86% of the sample is covered by the Spitzer observations 
down to a limiting flux of ~ 83 /i Jy. The detected galaxies 
yield a median ratio of (SFRir/SFRuv) of ~ 3.6, indica- 
tive of a sustantial amount of obscured star formation. 
Three of the interacting galaxies in the first redshift bin 
had anomalously high SFRuv+iR 41, 18, and 15 Mq 
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yr~^). Two of these turned to have infrared spectra con- 
sistent with an AGN and were removed before computing 
the IR-based SF properties shown in Fig. [T3]to Fig. 

The average UV-based SFR (based on 3698 galaxies) 
and UV+IR-based SFR (based on only the 876 galaxies 
with 24um detections) are plotted in the top 2 panels 
of Fig. [T3] for three groups of intermediate mass galax- 
ies: mergers, non-Interacting E-Sd, and non-Interacting 
Irrl. The corresponding plot for the high mass sample is 
in Fig. [ini It can be seen (Fig. [H] and Fig. [T5)) that over 
z ~ 0.24-0.80, the average UV-based and UV-|-IR-based 
SFR of mergers (in the phase where they are recognizable 
as mergers) are only modestly enhanced^ at best by a factor 
of a few, compared to the non-interacting galaxies. This 
result applies to both high mass (M > 2.5 x 10^° Mq) 
galaxies and intermediate mass (M > 1.0 x 10^ Mq) blue 
cloud galaxies. A similar result is also found by Robaina 
et al. (in preparation) in high mass systems. This mod- 
est enhancement is consistent with the recent statistical 
study of di Matteo et al. (2007), who find from numerical 
simulations of several hundred galaxy collisions that the 
maximum SFR in galaxy mergers is typically only a factor 
of 2-3 larger than that of corresponding non-interacting 
galaxies. Their results suggest that the results of some 
early simulations (e.g., Mihos & Hernquist 1996; Hern- 
quist & Mihos 1995), where mergers converted 50 to 80 
per cent of their original gas mass into stars, may not rep- 
resent the typical situation at z < 1. 

In order to further test the robustness of our result, 
we used the stacking procedure described in Zheng et 
al. (2006) to get a more representative measure of the IR- 
based SFR for the following three groups of intermediate 
mass systems: mergers, non-interacting E-Sd, and non- 
interacting Irrl galaxies. For every group, the individual 
galaxies were cross-correlated with the Spitzer 24 /im cat- 
alog in order to identify detected and undetected objects. 
Then the PSF-removed 24 /xm images for the undetected 
objects were stacked, and a mean flux was derived from the 
average/median stacked image. An average 24um luminos- 
ity was determined from the individually-detected fluxes 
and individually-undetected fluxes estimated by stacking. 
The 3215 intermediate mass galaxies in the Spitzer field 
were used in this process, giving a more representative 
24um luminosity than the mere 878 galaxies with detec- 
tions. A final uncertainty can be obtained by combining 
background error and bootstrap error in quadrature. The 
IR-based SFR was estimated from the 24um luminosity 
using the procedure described above, and combined with 
the UV-based SFR to estimate the total SFR. The aver- 
age UV+IR-stacked SFR is plotted in the bottom panel 
of Fig. [13] and Fig. [151 again, only a modest enhance- 
ment is seen in the average SFR of mergers (in the phase 
where they are recognizable as mergers) , compared to non- 
interacting galaxies. 

4.7. The contribution of interacting galaxies to the 
cosmic SFR density over the last 7 Gyr 

Over the last 8 Gyr since z ^ 1, the cosmic SFR density 
is observed to decline by a factor of 4 to 10 (e.g., Lilly et 
al. 1996; Ellis et al 1996; Flores et al. 1999; Haarsma 
et al. 2000; Hopkins 2004; Perez-Gonzalez et al. 2005; 
Le Floc'h et al. 2005). Earlier GEMS studies by Wolf 



et al. (2005) and BeU et al. (2005) over a 0.6 Gyr inter- 
val (z - 0.65-0.75 or Tb - 6.2-6.8 Gyr) showed that the 
UV and IR luminosity density over this interval are dom- 
inated by non-interacting galaxies. Here, we extend the 
earlier GEMS studies to cover a six-fold larger time inter- 
val of 4 Gyr {z ^ 0.24-0.80 or Tb ~ 3-7 Gyr), and set 
quantitative limits on the contribution of merging systems 
to the UV-based and UV+IR-based SFR density. We use 
the sample S2 of ^^3698 intermediate mass (M > 1.0 x 10^ 
Mq) galaxies. Our study also complements the IR-based 
studies by Hammer et al. (2005; 195 galaxies at z > 0.4 ) 
and Melbourne et al. (2005; ^ 800 galaxies) in terms of 
sample size or/and SFR indicators. 

Fig. [16] shows the SFR density for intermediate 
mass mergers, non-interacting E-Sd galaxies, and non- 
interacting Irrl galaxies over z ^ 0.24-0.80. The top panel 
shows the UV-based SFR density from the full sample. 
The middle panel show the UV-|-IR-based SFR density 
from the 878 galaxies with individual 24um detections. 
Finally, the bottom panel shows the UV-|-IR-stacked SFR 
density determined via the stacking of 3215 galaxies with 
Spitzer coverage, as outlined in § 14.61 In all three pan- 
els, one finds that interacting galaxies only account for a 
small fraction (< 30%) of the cosmic SFR density over 
z ~ 0.24~0.80, corresponding to lookback times of 3-7 
Gyr (Fig. [75)) . The same results hold for the sample of 
high mass (M > 2.5 x 10^*^ Mq) galaxies, as illustrated in 
Fig. [13 

Thus, our results suggest that the behavior of the cosmic 
SFR density over z ~ 0.24-0.80 is predominantly shaped 
by non- interacting galaxies. A similar result is found by 
Lotz et al. (2008). Our result is a direct consequence of 
the fact that the merger fraction / (Table [21 ii l4.3p . as well 
as the enhancement in the average SFR from interactions 
(§ 14. 6p . are both modest. Our results agree remarkably 
well with models for the self-regulated growth of super- 
massive black holes in mergers involving gas-rich galaxies 
(Hopkins et al. 2006). These models predict that galaxy 
mergers contribute only ~ 20% of the SFR density out to 
z ^ 1, and even out to z 2. 

It is legitimate to ask whether the results hold despite 
the uncertainties in identifying mergers. We first note that 
based on the tests of § 14.11 we have already included a 
large fractional error term on / to account for the bino- 
mial standard deviation, the dispersion between classifiers, 
and the effect of moderate bandpass shifting, and surface 
brightness dimming. Therefore, the results presented here 
already take into account at least some of these sources of 
uncertainties. 

Another source of uncertainty might be that some of 
the galaxies, which we have classified as non-interacting 
Irrl under the assumption that their small-scale asym- 
metries are likely caused by SF rather than interactions 
(§ 13. 3[) . may be borderline cases of mergers or interact- 
ing systems. However, it can be seen from Fig. 1161 that 
even if we were to add the SFR density of all the non- 
interacting Irrl galaxies to that of the mergers, the sum 
would still be significantly lower than the contribution of 
non-interacting E to Sd galaxies. Thus, the results would 
be largely unchanged. 

Another test is to repeat the analyses using the CAS 
merger criterion {A > 0.35 and ^ > S") to identify merg- 
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ers. The limited recovery rate (50% to 73%) and large level 
of contamination impacting the CAS criterion (ij |3.4p make 
it more difficult to interpret the SF properties of systems 
identified as mergers or non-interacting with CAS. For 
both the intermediate mass sample (Fig. [T8| and the high 
mass sample (Fig. [11]), the average SFR of CAS mergers is 
only modestly enhanced compared to CAS non-interacting 
galaxies, in agreement with the results from § 14.61 

Furthermore, for the intermediate mass sample. Fig. 1201 
shows that CAS mergers contribute only 16% to 33% of 
the UV SFR density and 22% to 38% of the VV+IR SFR 
density. While the upper limits of these values are slightly 
higher than those based on the visual types, it is nonethe- 
less reassuring that CAS non-interacting galaxies domi- 
nate the SFR density. Similar results hold for the sample 
of high mass {M > 2.5 x 10^° Mq) galaxies, as illustrated 
in Fig. (211 

For intermediate mass {M > 1.0 x 10^ -^0) galaxies, 
we find that the cosmic SFR density declines by a fac- 
tor of - 3 from z - 0.80 to 0.24 (lookback time ~ 7 to 
3 Gyr). Since non-interacting galaxies dominate the cos- 
mic SFR density in every redshift bin, it follows that this 
decline is largely the result of a shutdown in the SF of 
non-interacting galaxies. The question of what drives this 
shutdown will be addressed in detail in a future paper, and 
is only considered briefiy here. One possibility is the de- 
pletion of the internal cold gas supply of galaxies by star 
formation, or the reduction in the accretion rate of gas 
from cosmological filaments. Future facilities like ALMA 
will be instrumental in exploring this issue further. An- 
other related possibility is that over time, most of the SFR 
is shifting to lower stellar masses. High mass systems are 
associated with a lower SSFR (Cowie et al. 1996; Brinch- 
mann et al. 2004; Brinchmann & Elhs 2000; Fontana et 
al. 2003; Bauer et al. 2005; Zheng et al. 2007; Fig. 1 of 
Noeske et al. 2007a), consistent with the idea that they 
have experienced the bulk of their stellar mass growth at 
earlier epochs {z > 1). In staged SF models (Noeske et al. 
2007b) the SF history of low mass systems is consistent 
with exponential SF models associated with a late onset 
and a long duration. 

5. SUMMARY AND CONCLUSIONS 

We have performed a comprehensive observational es- 
timate of the galaxy merger fraction over z ~ 0.24-0.80 
(lookback times of 3-7 Gyr), and explored the impact of 
mergers on the star formation of galaxies over this inter- 
val. Our study is based on HST ACS, COMBO-17, and 
Spitzer 24 /im data from the GEMS survey. We use a large 
sample of - 3600 (M > 1 x 10^ Mq) galaxies and 790 
high mass (M > 2.5 x 10^" Mq) galaxies (§[!]). We pri- 
marily identify mergers using a visual classification system, 
which is based on visual morphologies, spectrophotometric 
redshifts, and stellar masses 13.21 to 13. 3p . and identi- 
fies systems that show evidence of having experienced a 
merger of mass ratio M1/M2 > 1/10 within the last vis- 
ibility timescale. While many earlier studies focused only 
on major mergers (defined as mergers with A/1/M2 > 
1/4), we also attempt to constrain the frequency of minor 
mergers (defined as mergers with 1/10 < M1/M2 < 1/4), 
since they dominate the merger rates in ACDM models. 
Below is a summary of our results: 



1. For the high mass {M > 2.5 x 10^" Mq) sample of 
^ 790 galaxies, which is complete on both the blue 
cloud and red sequence, we find the following. The 
fraction / of visually-classified systems that show 
evidence of a recent merger of mass ratio > 1/10, 
does not show strong evolution over lookback times 
of 3-7 Gyr, and ranges from 9% ± 5% at 2 - 0.24- 
0.34, to 8% ± 2% at z - 0.60-0.80 (Tabled Figg]). 
These mergers are further subdivided into three cat- 
egories: clear 'major merger', clear 'minor merger', 
and ambiguous 'major or minor merger' cases. The 
first two classes are used to set lower limits on the 
major and minor merger fraction. The lower limit 
on the major merger fraction, determined in this 
way, ranges from 1.1% to 3.5% over z ~ 0.24-0.80 
(Table[T]). The corresponding lower limit on the mi- 
nor merger fraction ranges from 3.6% to 7.5%. This 
is the first, albeit approximate, empirical estimate 
of the frequency of minor mergers over the last 7 
Gyr. 

For an assumed visibility timescale of ~ 0.5 Gyr, 
it follows that over Tb ~ 3-7 Gyr, - 68% of high 
mass systems have undergone a merger of mass ra- 
tio > 1/10, with ~ 16%, 45%, and 7% of these cor- 
responding respectively to major, minor, and am- 
biguous 'major or minor' mergers. The correspond- 
ing merger rate i? is a few xlO~^ galaxies Gyr~^ 
Mpc-3. 

2. At intermediate masses (M > 1 x 10^ Mq), we are 
only complete in mass for the blue cloud. Among 
- 2840 blue cloud galaxies of mass M > 1.0 x 10^ 
Mq, the fraction of visually-classified systems that 
show evidence of a recent merger of mass ratio 
> 1/10, ranges from 7% ± 2% to 15% ± 5% over 
z - 0.24-0.80 (Tabled]). 

For an assumed visibility time of ^ 0.5 Gyr, we es- 
timate that on average, over Tb ~ 3-7 Gyr, 84% of 
intermediate mass blue cloud galaxies have under- 
gone a merger of mass ratio > 1/10, with ^ 5%, 
22%, and 57% corresponding respectively to ma- 
jor, minor, and ambiguous 'major or minor' merg- 
ers. The corresponding merger rate R ranges from 
8 X 10"^ to 1 X 10-3 galaxies Gyr"! Mpc'^ . 

3. We compare our visual mergers to those identi- 
fied using the widely used CAS merger criterion 
{A > 0.35 and ^ > S'), based on CAS asymme- 
try A and dumpiness S parameters (§ 14. 2p . The 
merger fraction based on the CAS merger criterion 
agrees within a factor of two with the visually based 
merger fraction for high mass (M > 2.5 x 10^° 
Mq) galaxies. However, for intermediate mass 
(M > 1 X 10® Mq) galaxies, CAS can overestimate 
the merger fraction at z > 0.5 by a factor ^ 3. In ef- 
fect, over z ~ 0.24-0.80, - 50% to 70% of the galax- 
ies visually classified as mergers satisfy the CAS cri- 
terion, but the latter also picks up a dominant num- 
ber of non-interacting dusty, star-forming galaxies 
(Fig. [7] and Fig. [5]). These non-interacting systems 
make up as much as ~ 45% to 80% of the systems 
picked up the CAS criterion. We thus conclude that 
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the traditional CAS merger criterion is ill-suited for 
use on HST V-h&nd images at z > 0.5 where the 
rest frame wavelength falls below A < 4000 A, par- 
ticularly in the case of intermediate mass galaxies 
with significant SF, gas, and dust. Modified CAS 
criteria in the near-UV based on morphological k- 
corrections (Taylor et al. 2007) might alleviate this 
problem. 

4. We compare our empirical merger fraction / and 
merger rate R for high mass (M > 2.5 x 10^° 
Mq) galaxies to predictions from different ACDM- 
based simulations of galaxy evolution, including the 
halo occupation distribution (HOD) models of Hop- 
kins et al. (2007); semi-analytic models (SAMs) of 
Somerville et al. (2008), Bower et al. (2006), and 
Khochfar & Silk (2006); and smoothed particle hy- 
drodynamics (SPH) cosmological simulations from 
Mailer et al. (2006) with a corrected stellar mass 
function (see §113]). To our knowledge, such exten- 
sive comparisons have not been attempted to date, 
and are long overdue. We find qualitative agree- 
ment between the observations and models, such 
that the (major+minor) merger fraction fFia. [T/]) 
and merger rate (Fig. [71|j from different models 
bracket the corresponding empirical estimates and 
show a factor of five dispersion. One can now an- 
ticipate that in the near future, improvements in 
both the observational estimates and model predic- 
tions will start to rule out certain merger scenarios 
and refine our understanding of the merger history 
of galaxies. 

5. We explore the impact of galaxy mergers on the SF 
activity of galaxies since z < 0.8. In the sample of 
~ 789 high mass (M > 2.5 x 10^" Mq) galaxies, 
as well as the sample of ~ 3600 intermediate mass 
(M > 1.0 X 10^ Mq) galaxies, we find that the av- 
erage SFR of visibly merging galaxies is only mod- 
estly enhanced compared to non-interacting galax- 
ies over z ~ 0.24-0.80 (Fig. \T^. This result is 
found for SFRs based on UV, UV-I-IR, as well as 
UV-|-stacked-IR data. This modest enhancement is 
consistent with the empirical results of Robaina et 
al. (in preparation), and the recent statistical study 
of di Matteo et al. (2007) based on numerical sim- 
ulations of several hundreds of galaxy collisions. 

6. Among both high mass and intermediate mass 
galaxies, our results of a modest merger fraction 
/ and a modest enhancement in the average SFR 
due to mergers, culminate in our finding that visibly 
merging systems only account for a small fraction 
(< 30%) of the cosmic SFR density over lookback 
times of ~ 3-7 Gyr (z ~ 0.24-0.80; Fig. \M and 
Fig. [77y . Our result complements that of Wolf et 
al. (2005) over a smaller lookback time interval of 
~ 6.2-6.8 Gyr. In effect, our result suggests that 
the behavior of the cosmic SFR density over the last 
7 Gyr is predominantly shaped by non-interacting 
galaxies, rather than mergers and interacting galax- 
ies. We suggest that our observed decline in the cos- 
mic SFR density by a factor of ~ 3 since z ~ 0.80 



is largely the result of a shutdown in the SF of rela- 
tively non-interacting galaxies. This shutdown may 
be driven by the depiction of the internal cold gas 
supply of galaxies, the reduction in the accretion 
rate of gas from cosmological filaments, and the 
transition of SF activity to lower mass systems. 
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Fig. 1. — The rest-frame U — V color is plotted versus the stellar mass over z 0.24-0.80 for the sample S2 of ~ 3698 galaxies with 
M > 1 X 10^ Mq. The four panels denote the four redshift bins, which span 1 Gyr each, and cumulatively cover the interval z ~ 0.24-0.80 
(Tback ~ 3—7 Gyr). N denotes the number of galaxies in each bin. The diagonal line marks the separation of the red sequence galaxies and 
the blue cloud galaxies at the average redshift Zave of the bin. The vertical lines marks the mass completeness limit (Borch et al. 2006) for the 
red sequence galaxies, while the blue cloud galaxies are complete well below this mass. For the mass range M > 1 X 10^ A/q, the blue cloud 
is complete in our sample out to z ~ 0.80, while the red sequence is incomplete in the higher redshift bins. Galaxies are coded according to 
their visual type (VT) in the F606W band: 'Mergers' (orange stars), 'Non-Interacting E+SO-|-Sa' (black diamonds), 'Non-Interacting Sb-Sc 
-I- Sd' (green diamonds), and 'Non-Interacting Irrl' (blue squares). 



Galaxy Interactions and Star Formation over 7 Gyr 



19 




4 2=0.42 5^ 2=0.43 6 2=0.44 
7 2=0.53 e 2=0.56 9 I 2=0.66 




Fig. 2. — This montage show examples of mergers, namely systems that show evidence of having experienced a merger of mass ratio 
> 1/10 within the last visibility timescale, as described in ^ 13.21 The evidence includes morphological distortions similar to those seen in 
simulations of mergers of mass ratio > 1/10, such as multiple nuclei (e.g., case 6, 8); components (e.g., case 12) connected by a bridge or 
common envelope; tidal tails and asymmetric features (e.g., cases 3, 4, 5, 9, 11); and warped disks (e.g., case 2). Systems classified as 'Int-1' 
mergers primarily represent advanced mergers, which appear as single systems in ACS images (e.g., cases 2, 3, 4, 5, 6, 7, 9, 10 11, and 12). 
Conversely, 'Int-2' mergers (e.g., case 1) primarily represent young mergers, which appear as very close pairs of overlapping galaxies in ACS 
images. The mergers can be divided into 3 groups: clear major mergers (cases 1, 6, 12), clear minor mergers (cases 2, 9), and ambiguous 
'major or minor mergers' (cases 3, 4, 5, 7, 8, 10, 11). 
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Fig. 3. — This montage shows examples of galaxies classified as non- Interacting E-to-Sd and non-interacting Irrl galaxies, according to the 
criteria in § 13.31 Within the broad class of non-interacting E-to-Sd, galaxies have Hubble types E, SO, Sa, Sb-Sc, and Sd, as shown in Fig. |6] 
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Fig. 4. — This figure compares the the merger fraction (/) based on visual classification by 3 classifiers (SJ, SM, KP), to merger fraction 
(/cAs) that would be obtained using the CAS criterion {A > 0.35 and A > S). The results are shown for both high mass {M > 2.5 X 10^" 
Mq; top panel) and intermediate mass {M > 1 X 10^ ^Q'l bottom panel) sample. The plotted error bar for /, at this stage, only includes the 
binomial term [/(l-/)/A'^]^/^, for each bin of size A^. The same trend is seen for all 3 classifiers and the maximum spread Sf/f in the 4 bins is 
~ 26%. /cAS agrees within a factor of two with the visually based / merger fraction for high mass galaxies. However, for intermediate mass 
galaxies, CAS can overestimate the merger fraction at z > 0.5 by a factor ~ 3, as it picks up a significant number of non-interacting dusty, 
star-forming galaxies (see 14.21 1. 
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Fig. 5. — This montage illustrates a test for bandpass shift and surface brightness dimming. It compares the bluer shallower GEMS 
F606W images (V band; pivot A ~ 5915 A) and deeper redder GOODS F850LP (z band; pivot A ~ 9103 A) images of typical interacting and 
non-interacting galaxies in the last redshift bin (z ~ 0.60 to 0.80), where bandpass shift and surface brightness dimming are expected to be 
most severe. In this redshift bin, the rest-frame wavelength traced by the GEMS images shift from optical to violet /near-UV (3700 A to 3290 
A). However, while the GOODS images have higher S/N, and trace redder older stars, they do not reveal dramatically different morphologies 
from from those in the GEMS F606W images. Cases 1, 5, 6, 7, and 8 are interacting systems. A statistical analysis is shown in Table [3] 
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Fig. 6. — For intermediate mass (M > 1 X 10^ Mq) galaxies, the distribution of Sersic index n from single-component Sersic fits is plotted 
for mergers, non-interacting Irrl, and non-interacting E-to-Sd systems. The latter class is further subdivided as E, SO, Sa, Sb-Sc, and Sd. 
The majority of systems visually classified as non-interacting Sb-Sc, Sd, and Irrl have n < 2.5, as expected for disk-dominated systems. Most 
of the systems visually typed as Sa have n < 4. Those galaxies typed as E and SO primarily have n > 3, but some have low n, reflecting the 
inherent difficulty in separating E, SO, and Sa at intermediate redshifts. 
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Fig. 7. — The CAS asymmetry A and dumpiness S parameters are plotted for intermediate mass (A/ > 1 X 10^ ^^0) galaxies in the 
four redshift bins of Fig. [T] using the same color coding. Galaxies satisfying the CAS criterion {A > 0.35 and A > S) lie in the upper left 
hand corner, bracketed by the A = S diagonal line and the A = 0.35 horizontal line. The CAS criterion captures a fair fraction of the galaxy 
mergers, but it also picks up "contaminants" in the form of non-interacting galaxies. This is further illustrated in Fig. \E\ 
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Fig. 8. — For galaxies with M > 1.0 X 10^ Mq, the top panel shows that the fraction i^cAS— merger of visually-classified mergers, which 
satisfy the CAS criterion is ~ 50% to 70% across the four redshift bins. The bottom panel shows the degree to which non-interacting 
galaxies contaminate the systems picked up by CAS. A^CAS represents the total number of galaxies, which satisfy the CAS criterion and are 
considered as "CAS mergers" across the four redshift bins. The fraction i^cAS-visual of different visual types among these "CAS mergers" is 
plotted on the y-axis. At 2 > 0.5, the vast majority (44% to 80%) of the systems considered as mergers by the CAS criterion turn out to be 
non-interacting [E-Sd and Irrl] systems. 
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Fig. 9. — The montage shows typical systems where the CAS criterion (A > 0.35 and A > S) fails (see ij 14.21 for details). Cases 1-3 are 
systems, which are visually classified as mergers, but are missed by the CAS criterion. They include systems with tidal debris (e.g., case 3) 
that may contribute less than 35% of the total light; systems with close double nuclei (e.g., case 2) where CAS might refine the center to 
be between the two nuclei, thereby leading to a low A < 0.35; and pairs of fairly symmetric galaxies whose members have similar redshifts 
within the spectrophotometric error, appear connected via weak tidal features, and have a stellar mass ratio M1/M2 > 1/10 (e.g., case 1 
where MI/M2 ~ 0.25). Conversely, cases 4-9 are systems, which are visually classified as non-interacting galaxies, but are picked by the CAS 
criterion. They include non-interacting, actively star- forming systems with small-scale asymmetries in the optical blue light (cases 4 and 6); 
systems where A is high due to the absence of a clear center (case 8) or due to the center being blocked by dust (case 4, 9); edge-on systems 
and compact systems, where the light profile is steep such that small centering inaccuracies lead to large A (case 9). 
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Fig. 10. — The observed merger fraction / in the high mass (M > 2.5 X 10^" Mq) sample is compared to other studies, noting the caveats 
outlined in ii 14.41 Shown here are the merger fraction based primarily on morphologically distorted galaxies (filled circles: this study; Lotz et 
al. 2008; Conselic e 20 03), and the close pair fraction (open squares: Le Fevre et al. 2000; Bell et al. 2006; Kartaltepe et al. 2007) as a function 
of redshift. See ij l4.4l for details. 



28 



Jogee, Miller, Penner, and the GEMS collaboration 



Massive (M/Mo>2.5x1 0'°) Galaxies [N = 789] 



c 
o 

o 
o 



0) 
0) 



OJO 




Data:.)K (Mojor+Minor) 

Q Q1 Models: (Mojor+Mlnor) (Mojor) 

Bower 06 (SAM) Khochfor 06 (SAM) Somerville 08 (SAM) 
Moller 06 (SPH) Hopkins 07 (HOD) 



0.0 



0.2 



0.4 

Redshift 



0.6 



0.8 



Fig. 11. — The empirical merger fraction / (orange stars) for mergers with mass ratio M1/M2 > 1/10 among high mass galaxies is compared 
to the fraction of (major+minor) mergers (solid lines; stellar mass ratio M1/M2 > 1/10) predicted by different ACDM-based simulations of 
galaxy evolution, including the halo occupation distribution (HOD) models of Hopkins et al. (2007); semi-analytic models (SAMs) of Somerville 
et al. (2008), Bower et al. (2006), and Khochfar &i Silk (2006); and smoothed particle hydrodynamics (SPH) cosmological simulations from 
Mailer et al. (2006) (see ^{tb^. 



Galaxy Interactions and Star Formation over 7 Gyr 



29 



0.0010 



Massive (M/Mo>2.5x1 0^°) Galaxies [N = 789] 



I 



O 
Q_ 



I 

<u 0,0001 
-I—' 

o 












_ . — ■ ' 





Data: (Major+Minor) 

Models: (Mojor + Mlnor) 

Bower 06 (SAM) Khochfor 06 (SAM) 
Moller 06 (SPH) Hopkins 07 (HOD) 



(Major) 

Somerville 08 (SAM) 



0.0 



0.2 



0.4 
Redshift 



0.6 



0.8 



Fig. 12. — As in Figure 11, but now comparing the empirical rate R (orange stars) of mergers with mass ratio M1/M2 > 1/10 among 
high mass galaxies to the rate of (major+minor) mergers (solid lines; stellar mass ratio MI/M2 > 1/10) predicted by different ACDM-based 
simulations of galaxy evolution. 
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13. — The UV-based star formation rate is plotted versus the stellar mass over z ~ 0.24-0.80 for the sample S2 of galaxies with 
X 10^ Mq. The four panels show the four redshift bins, which span 1 Gyr each, and cumulatively cover the interval z ~ 0.24-0.80 
~ 3—7 Gyr). N denotes the number of galaxies plotted in each bin. Galaxies are coded as in Fig.[T] with merging systems denoted by 
stars. The average SFR and total SFR density in both the UV and the IR, are further illustrated in Fig. I14l and Fig. 1161 
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Fig. 14. — For the sample S2 of galaxies with M > 1 X 10^ Mq, the average SFR of merging systems, non-interacting E-Sd galaxies, and 
non-interacting Irrl galaxies are compared over z ^ 0.24-0.80. A'^ denotes the number of galaxies used. The average UV-based SFR (top 
panel; based on 3698 galaxies), average UV-I-IR- based SFR (middle panel; based on only the 876 galaxies with 24um detections), and average 
UV-|-IR-stacked SFR (based on 3215 galaxies with 24um coverage) are shown. In all there cases, the average SFR of visibly merging systems 
is only modestly enhanced compared to non-interacting galaxies over z ~ 0.24-0.80 (lookback time ~ 3-7 Gyr). See 14.61 for details. 
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Fig. 15. — As in Fig. 1141 but for the sample SI of high mass (M > 2.5 X 10^" Mq) galaxies. Only data points with at least 3 galaxies are 
shown. Again, the average SFR of visibly merging galaxies is only modestly enhanced compared to non-interacting galaxies over z ~ 0.24-0.80 
(lookback time ~ 3—7 Gyr). 
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Fig. 16. — For the sample S2 of galaxies with M > 1 X 10® Mq, the SFR density of merging systems, non-interacting E-Sd galaxies, 
and non-interacting Irr galaxies are compared over z ^ 0.24—0.80. Results based on UV (top panel), UV+IR (middle panel), as well as 
UV-|-stacked-IR data (bottom panel), are shown in the top, middle, and bottom panels. In all bins, visibly merging systems only contribute 
a small fraction (typically below 30%) of the total SFR density. In effect, the behavior of the cosmic SFR density over the last 7 Gyr is 
predominantly shaped by non-interacting E-Sd galaxies rather than visibly merging galaxies. 
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Fig. 17. — As in Fig. 1161 but for the sample of high mass (Af > 2.5 X 10^" Mq) galaxies. Only data points with at least 3 galaxies are 
shown. The same conclusion holds: the cosmic SFR density over the last 7 Gyr is predominantly shaped by non-interacting E-Sd galaxies 
rather than visibly merging galaxies. 
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Fig. 18. — Same as in Fig. 1141 but using the CAS criterion {A > 0.35 and A > S) to identify interacting systems in the sample of 
intermediate mass (M > 1.0 X 10^ Mq) galaxies. The average SFR of ' C AS-interacting' galaxies is only modestly enhanced compared to 
'CAS non-interacting' galaxies, in agreement with the results from ^ 14.61 
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Fig. 19. — As in Fig. [Tsl but for the sample of high mass (M > 2.5 X 10^0 Mq) galaxies. The same conclusion holds: the average SFR of 
'CAS-interacting' galaxies is only modestly enhanced compared to 'CAS non-interacting' galaxies. 
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Fig. 20. — Same as in Fig. 1161 but using the CAS merger criterion {A > 0.35 and A > S) to identify interacting galaxies in tlie sample of 
intermediate mass (A/ > 1.0 X 10^ Mq) galaxies. The 'CAS-interacting' galaxies contribute only 16% to 33% of the UV SFR density and 22% 
to 38% of the UV+IR SFR, density. While the upper limits of these values are slightly higher than those based on the visual typos (Fig. I16|l. 
the 'CAS non-interacting' galaxies' clearly dominate the SFR, density. 
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As in Fig. 1201 but for the sample of high mass (M > 2.5 X 10^" Mq) galaxies. The same conclusion holds. 
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Table 1 

YisvxL AlERCER Fraction for Hicii Mass (M 2.5 x 10"^" M . ) Sample SI [X-7891 



(1) 


Redshift bin 


1 


2 


3 


4 


(2) 


Redshift range 


0.24-0.34 


0.34-0.47 


0.47-0.62 


0.62-0.80 


(3) 


Lookback time [Gyr] 


3.0-4.0 


4.0-5.0 


5.0-.6-0 


6.0-7.0 


(4) 


Arest in F606W [A] 


4470-4414 


4414-4023 


4023-3651 


3651-3286 


(5) 


Total no of galaxies 


46 


84 


213 


446 


(6) 


Fraction of mergers 


0.087±0.047 


0.083±0.037 


0.089±0.030 


0.0807±0.025 


(6a) 


Lower limit on major merger fraction 


0.022±0.021 


0.035±0.022 


0.014±0.009 


0.011±0.006 


(6b) 


Lower limit on minor merger fraction 


0.065±0.040 


0.036±0.022 


0.075±0.027 


0.049±0.016 


(6c) 


Fraction of ambiguous minor /major mergers 


0.00 


0.012±0.012 


0.00 


0.020±0.008 


(7) 


Fraction of non-interacting E-Sd 


0.913± 0.241 


0.869±0.229 


0.878±0.229 


0.785±0.205 


(8) 


Fraction of non-interacting Irrl 


0.000 


0.024±0.018 


0.009±0.007 


0.025±0.010 


(9) 


Fraction of compact 


0.000 


0.024±0.018 


0.023±0.012 


0.11±0.032 



Note. — Rows are : (1) Redshift bin. (2) Range in redshift covered by the bin; (3) Range in lookback time covered by the bin; (4) Range 
in rest-frame wavelength traced by the F606W filter over the bin, assuming a pivot wavelength of 5915A; (5) Total number of high mass 
galaxies per bin; (6) Fraction of systems with evidence of a recent merger of mass ratio > 1/10. These include both major {M1/M2 > 1/4) 
and minor (1/10< MI/M2 < 1/4) mergers; (6a) Lower limit on the fraction of galaxies undergoing major mergers; (6b) Lower limit on the 
fraction of galaxies undergoing minor mergers; (6c) Remaining fraction of galaxies that could be either major or minor mergers; (7-9) Fraction 
of non-interacting E-Sd, non-interacting Irrl, and compact systems. 
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Table 2 

Visual Merger Fraction for M > 1.0 x 10^ Mq Sample S2 [N=3698] 



(1) 
(2) 
(3) 
(4) 



Redshift bin 
Redshift range 
Lookback time [Gyr] 
Arest in F606W [A] 



0.24-0.34 

3.0-4.0 
4470-4414 



0.34-0.47 
4.0-5.0 
4414-4023 



0.47-0.62 

5.0 .6-0 
4023-3651 



0.62-0.80 

6.0-7.0 
3651-3286 



All [N=3698] 

(5) Total no of galaxies 235 480 1117 1866 

(6) Fraction of mergers 0.111±0.035 0.090±0.027 0.074±0.021 0.062±0.017 
(6a) Lower limit on major merger fraction 0.013±0.008 0.006±0.004 0.004±0.002 0.009±0.003 
(6b) Lower limit on minor merger fraction 0.038±0.016 0.021±0.008 0.021±0.007 0.017±0.005 
(6c) Fraction of ambiguous minor/major mergers 0.060±0.022 0.062±0.020 0.048±0.014 0.036±0.010 

(7) Fraction of non-interacting E-Sd 0.850±0.220 0.846±0.220 0.796±0.207 0.793±0.206 

(8) Fraction of non-interacting Irrl 0.064±0.023 0.052±0.017 0.108±0.030 0.064±0.018 

(9) Fraction of compact 0.000 0.012±0.006 0.021±0.007 0.080±0.022 



Blue Cloud [N=2844] 

(10) Total no of galaxies 154 332 876 1482 

(11) Fraction of mergers 0.149±0.048 0.114±0.034 0.088±0.025 0.069±0.019 
(11a) Lower limit on major merger fraction 0.013±0.009 0.00 0.005±0.003 0.008±0.003 
(lib) Lower limit on minor merger fraction 0.046±0.020 0.024±0.010 0.023±0.008 0.016±0.005 
(He) Fraction of ambiguous minor/major mergers 0.091±0.033 0.090±0.028 0.060±0.018 0.046±0.013 

(12) Fraction of non-interacting E-Sd 0.753± 0.199 0.801±0.209 0.753±0.196 0.784±0.204 

(13) Fraction of non-interacting Irrl 0.097±0.035 0.075±0.024 0.136±0.037 0.080±0.022 

(14) Fraction of compact 0.000 0.009ib0.006 0.023ib0.008 0.067ib0.018 



Note. — Rows arc : (1) to (9): As in Table 1, but for intermediate mass (M > 1.0 X 10^ ^^q) galaxies. However, note that the intermediate 
mass sample is incomplete for the red sequence. (10) to (14); Ditto, but for the blue cloud, where the sample is complete. 
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Table 3 

Merger fraction in GEMS F606W {V) and GOODS F850LP (z) [Ar=855] 







GEMS V 


GOODS z 


GOODS z 


GOODS z 


GOODS z 






Average 


Average 


S3 


SM 


KP 


(1) 


Fraction / of merging galaxies 


0.046±0.007 


0.049±0.007 


0.051±0.007 


0.057±0.008 


0.038±0.006 


(2) 


Ratio of / in GOODS z to GEMS V 




1.06±0.22 


1.10 ±0.23 


1.20±0.25 


0.83 ±0.18 



Note. — As a test for bandpass shift and surface brightness dimming, the table shows a comparison of the frax;tion of visual mergers based 
on GEMS F606W images and deeper redder GOODS F850LP images. The sample consists of the 855 intermediate mass (M > 1.0 x 10^ Mq) 
galaxies at z ~0.24 to 080, which are common to both GEMS F606W and GOODS F850LP surveys. Columns are: (2) Fraction / of mergers 
in GEMS F606W. The error bar only includes the binomial term (3) Average fraction / of mergers in GOODS F850LP based 

on results by classifiers (SJ,SM,KP), shown in columns 4-6. 



